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ABSTRACT 


This  thesis  introduces  a  planning  system  that  decides  which  missiles  to  load  on 
which  deploying  ships  to  maximize  their  campaign  effectiveness  across  all  anticipated 
operational  theater  war  plans.  Currently,  operational  planners  manually  identify  missile 
loadouts  and  ship  assignments  with  little  to  no  metric  to  identify  if  a  better  plan  exists. 
The  load  plans  must  be  robust  with  respect  to  a  range  of  potential  missions  and  conflicts 
(we  do  not  get  to  choose  the  war  we  must  fight)  and  must  provide  adequate  defensive 
ability  for  each  ship  and  offensive  ability  of  groups  of  ships.  We  consider  about  30-40 
customer  combatants  with  Mark  41  (MK41),  Mark  57  (MK57),  or  any  other  Vertical 
Launching  System  (VLS)  variant  to  appear,  including  ships  home-ported  in  the  U.S.  as 
well  as  forward-deployed  units.  There  are  nine  types  of  missiles,  limits  on  the  numbers  of 
each  missile  type  available,  the  amount  of  swapping  of  missiles  between  ships, 
compatibility  of  certain  types  with  certain  VLS  modules  on  certain  ships  and  perhaps 
concern  to  spread  certain  missiles  equitably  among  ships.  The  Excel-based  planning 
system  is  fast,  easy  to  understand  and  use,  and  operates  on  Navy  Marine  Corps  Intranet 
(NMCI)  computers. 
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EXECUTIVE  SUMMARY 


The  Mark  41  (MK  41)  Vertical  Launching  System  (VLS)  is  capable  of  holding 
offensive  and  defensive  weapons  with  numerous  combinations  of  load  outs.  Choosing  the 
right  load  out  for  a  given  mission  is  a  problem  for  operational  planners,  especially  since 
we  cannot  predict  a  priori  which  missions  will  be  required.  We  need  a  mission-robust 
load  out  plan — one  that  will  perform  well  across  a  broad  range  of  defenses  and  attacks. 
Commander  Seventh  Fleet  (C7F)  operational  planners  must  consider  a  number  of 
warplans  each  having  its  own  set  of  missions  and  a  number  of  ships  within  the  Area  of 
Responsibility  (AOR). 

This  problem  is  currently  being  addressed  by  operational  planners  using  basic 
Microsoft  Excel  (2010)  spreadsheets  to  provide  load  out  recommendations.  This  is  labor- 
intensive  and  provides  no  measurement  of  how  much  the  plans  might  be  improved. 

Wiederholt  (2015)  developed  and  implemented  a  mixed  integer  linear  program 
optimization  tool  known  as  the  Vertical  Launching  System  Loadout  Planner  (VLP)  in 
support  of  Seventh  Fleet  operational  planners  in  the  Western  Pacific.  This  thesis 
continues  his  work  and  matures  the  model.  The  product,  the  Missile  Loadout  Planner 
(MLP),  is  a  heuristic  solver  written  in  Visual  Basic  with  an  Excel  front  end  that  can  run 
on  Navy  and  Marine  Corps  Intranet  (NMCI)  and  Secret  Internet  Protocol  Router  Network 
(SIPRNet).  This  product  can  assist  operational  planners  in  the  daunting  task  of  assigning 
missile  loadouts  to  ships  in  preparation  for  deployment. 

MLP  is  tested  against  the  same  fictitious  warplans  as  the  Wiederholt  (2015) 
thesis.  There  are  two  warplans  and  52  missions  across  two  deployment  cycles.  The 
scenarios  use  23  VLS-capable  warships  and  nine  types  of  missiles.  The  planner  mimics 
the  mission  restrictions  of  VLP  and  suggests  a  missile  loadout  and  mission  assignment 
for  each  ship.  In  each  scenario,  the  ship  loadout,  mission  assignment,  and  missions 
unable  to  be  covered  due  to  shortages  are  provided.  The  scenarios  move  from  most 
restrictive  to  least  restrictive  with  regards  to  missile  loadout.  We  start  with  fixed  missile 
loadouts  and  only  recommend  mission  assignments  based  on  these  current  loadouts,  and 
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we  end  with  complete  flexibility  for  the  planning  tool  to  also  recommend  missile  loadout 
while  limiting  missile  movements  to  best  cover  each  mission  requirement.  This  tool  can 
assist  planners  by  providing  recommendations  in  a  fraction  of  their  current  manual 
planning  time.  Additionally,  this  tool  is  compatible  with  software  available  on  Navy 
network  systems. 
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INTRODUCTION 


The  Maritime  Operational  Planner  (MOP)  Office  of  Naval  Research  (ONR) 
Research  Program  at  the  Naval  Postgraduate  School  (NPS)  aims  to  produce,  develop, 
evaluate,  and  deliver  complete  operational  prototypes  of  various  optimization-based 
decision  support  tools  to  support  maritime  planning  staff  efforts.  This  thesis  is  part  of  the 
progression  of  NPS’s  larger  operations  research  program  to  provide  operational  planning 
decision  aids  for  the  maritime  operational  planner. 

Wiederholt  (2015)  began  the  process  here  by  investigating  and  implementing  an 
optimization  model  to  choose  the  best  missile  loadout  for  deploying  and  forward- 
deployed  ships  in  the  Commander  Seventh  Fleet  (C7F)  Area  of  Responsibility  (AOR). 
He  employed  proprietary  optimization  and  modeling  software  that  is  not  available  on 
Navy  and  Marine  Corps  Intranet  (NMCI)  and  in  some  facilities  served  by  Secret  Internet 
Protocol  Router  Networks  (SIPRNet) 

The  goal  of  this  thesis  is  to  continue  his  work  and  mature  the  model  to  achieve  a 
deployable  operational  planning  tool  that  all  fleets  and  naval  assets  can  use.  The  Missile 
Loadout  Planner  (MLP)  is  a  fast  heuristic  solver  written  in  Visual  Basic  for  Applications 
(Microsoft  2010)  with  a  Microsoft  Excel  2010  front  end  that  can  run  on  NMCI  and 
SIPRNet. 

A.  BACKGROUND 

This  section  will  provide  an  overview  of  the  problem.  It  will  cover  a  more 
detailed  description  of  the  problem  and  the  status  of  the  solution  as  a  result  of 
Wiederholt’ s  (2015)  work.  This  section  will  also  lay  out  some  of  the  foundation  of  the 
problem  including  detail  about  U.S.  Naval  missions,  ship  missile  systems,  and  ship 
missile  descriptions. 

1.  Problem  Description 

Ships  exist  in  two  main  states:  in  port  and  deployed.  Operational  planners  must 
choose  the  right  loadout  in  port  for  a  set  of  missions  to  be  executed  while  deployed.  The 
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MK  41  Vertical  Launching  System  (VLS)  is  capable  of  holding  offensive  and  defensive 
weapons  in  numerous  combinations  of  loadouts  (British  Aerospace  Engineering  [BAE] 
Systems  2017).  This  flexibility  can  only  be  used  to  our  advantage  if  planners  predict  a 
priori  which  missions  will  be  required.  Otherwise,  a  ship  might  enter  a  fight  without  the 
right  munitions.  This  problem  is  currently  addressed  using  basic  Microsoft  Excel  2010 
spreadsheets  to  evaluate  and  illustrate  loadout  recommendations.  This  method  is  labor- 
intensive  and  provides  no  metric  as  to  how  or  whether  the  plans  might  improve  given  a 
different  loadout.  Often,  a  ship’s  loadout  is  based  solely  on  which  missiles  are  coming 
out  of  maintenance  and  fleet  commanders  have  little  influence.  Therefore,  the  United 
States  Navy  (USN)  needs  a  mission-robust  loadout  plan  that  will  perform  well  across  a 
broad  range  of  anticipated  defensive  and  offensive  actions.  The  solution  to  this  problem 
is  important  because  it  is  vital  to  mission  success. 

The  missile  assignment  problem  is  further  complicated  by  other  restrictions.  For 
example,  the  cost  incurred  to  move  missiles  from  inventory  to  a  ship  or  to  switch  missiles 
between  ships  is  a  possible  deterrent  for  assigning  an  optimal  loadout.  Can  sensible 
reasons  be  found  to  govern  the  movement  of  missiles?  Second,  some  VLS-capable  ships 
cannot  carry  some  missile  types.  While  all  VLS  ships  can  carry  the  right  load  for  an 
escort  mission,  not  all  are  able  to  carry  Ballistic  Missile  Defense  (BMD)  loads.  In  some 
cases,  there  may  be  a  substitute  missile  for  a  mission  available  when  the  preferred  one  is 
not.  These  missile  substitutions  are  continually  increasing  as  newer  generations  of 
missiles  are  developed;  however,  the  VLS  itself  could  become  a  constraint.  Next- 
generation  missiles  have  the  potential  to  outgrow  the  size  limitations  of  the  current  MK 
41  VLS.  Last,  not  all  ships  are  in  port  at  one  time.  Some  ships  are  forward-deployed  at  all 
times.  These  ships  operate  in  conjunction  with  deploying  ships  from  ports  in  the  U.S. 
While  these  ships  are  deployed,  planners  are  already  preparing  the  loadout  for  the  next 
set  of  ships  that  will  act  as  a  relief  for  those  currently  in  theater. 

Wiederholt  (2015)  developed  and  implemented  a  mixed  integer  linear  program 
optimization  tool  known  as  the  VLS  Load  Planner  (VLP).  This  decision- support  tool 
determines  the  best  missile  loadout  that  can  be  achieved.  Although  his  research  supported 
Seventh  fleet  operational  planners  in  the  Western  Pacific  to  include  Forward  Deployed 
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Naval  Forces  (FDNF)  and  deploying  forces  to  the  area,  the  tool  can  be  applied  to  all  U.S. 
Fleet  AORs.  His  work  demonstrates  the  need  for  a  better  tool  than  current  methods  by 
showing  the  number  of  additional  missions  we  can  potentially  fulfill.  While  the 
Wiederholt  solution  is  optimal,  it  is  computed  using  software  (GAMS  2015)  not  available 
on  NMCI  or  all  SIPRNet  systems. 

2.  Department  of  Defense  Naval  Missions 

The  Department  of  Defense  (DOD)  issues  guidance  to  specify  definitions  of 
military  and  associated  terms.  This  thesis  will  use  the  following  definitions  to  describe 
the  missions  throughout: 

•  Air  defense  (AD):  Defensive  measures  designed  to  destroy  attacking 
enemy  aircraft  or  missiles  in  the  atmosphere,  or  to  nullify  or  reduce  the 
effectiveness  of  such  attack  (Joint  Chiefs  of  Staff  2016). 

•  Anti-submarine  warfare  (ASW):  Operations  conducted  with  the 
intention  of  denying  the  enemy  the  effective  use  of  submarines  (Joint 
Chiefs  of  Staff  2016). 

•  Protective  escort:  Defensive  posture  to  deter,  detect,  prevent,  and  defend 
against  attacks  to  high  value  units  (HVU)  (Department  of  the  Navy  [DoN] 
2014). 

•  Strike:  An  attack  to  damage  or  destroy  an  objective  or  a  capability  (Joint 
Chiefs  of  Staff  2016). 

•  Surface  warfare  (SUW):  That  portion  of  maritime  warfare  in  which 
operations  are  conducted  to  destroy  or  neutralize  enemy  naval  surface 
forces  and  merchant  vessels  (Joint  Chiefs  of  Staff  2016). 

•  Theater  Ballistic  Missile  Defense  (TBMD):  Primarily  defensive 
deployment  of  AD  missiles  to  protect  allies  in  a  specific  region,  or  theater 
(. Encyclopedia  Britannica  2016). 

3.  Vertical  Launching  System  Missiles 

The  following  missiles  and  descriptions  are  representative  of  those  used  in  the 
original  Wiederholt  formulation  and  include  updates  since  that  time.  The  three  letter 
missile  designation  describes  the  launch  environment,  mission  symbol,  and  type  of 
missile.  The  first  character,  “R,”  designates  a  ship-launched  missile;  the  second  character, 
“G,”  “I,”  or  “U”  indicates  a  surface  attack  missile,  aerial  intercept  missile,  or  underwater 
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attack  missile  respectively;  the  third  character,  “M”  indicates  a  guided  missile 
(Federation  of  American  Scientists  [FAS]  1998). 

•  RGM-109  Tomahawk  Land  Attack  Missile  (TLAM):  Surface  ship- 
launched  long  range,  subsonic  cruise  missile  used  for  strike  warfare 
missions  (FAS  2016). 

•  RIM-66  Standard  Missile  2  (SM2)  Medium  Range  (MR):  Includes 
Block  III,  IIIA,  and  IIIB  variants.  Primary  ship-launched  surface-to-air 
AD  and  ship  self-defense  missile  guided  by  semi-active  radar  or  infrared 
sensor  for  terminal  guidance  (United  States  Navy2016b). 

•  RIM-67  SM2  Extended  Range  (ER):  Now  known  as  the  RIM- 156  SM2 
Block  IV  variant.  Ship-launched  defensive  interceptor  to  provide  extended 
range  and  improved  high-altitude  air-defense  capability  (Global  Security 
2011). 

•  RIM-161  SM3:  Ship-launched  defensive  interceptor  to  counter  short  to 
medium-range  ballistic  missile  threats  (Raytheon  2017a). 

•  RIM-174  SM6  Extended  Range  (ER):  Ship-launched  over-the-horizon 
AD  and  SUW  missile  with  an  active  seeker  (United  States  Navy2016b). 

•  RIM-162  Evolved  Sea  Sparrow  Missile  (ESSM):  Ship-launched 
medium-range,  semi-active  homing  missile  for  ship  self-defense  against 
enemy  air  and  surface  threats  (United  States  Navy 2017a). 

•  RUM-139  Anti-submarine  rocket  (ASROC):  Surface  ship-launched 
anti-submarine  torpedo  (United  States  Navy 2013). 

4.  MK  41  Vertical  Launching  System  (VLS) 

The  MK  41  VLS  consists  of  multiple  eight-cell  launcher  modules  with  the  ability 
to  launch  SM  variants,  TLAM,  ASROC,  and  ESSM  in  support  of  multiple  missions. 
Figure  1  shows  an  example  of  one  such  module.  This  module  is  part  of  the 
Aegis  Weapons  System  and  is  installed  onboard  Arleigh-Burke  Class  Guided 
Missile  Destroyers  (DDG  51),  shown  in  Figure  2,  and  Ticonderoga  Class 
Guided  Missile  Cruisers  (CG  47),  shown  in  Figure  3,  (United  States  Navy2017b). 
The  modules  are  housed  below  the  deck  of  the  warship  to  protect  the  missiles  and 
make  the  system  more  survivable.  Each  module  has  the  ability  to  prepare  two 
missiles  simultaneously.  This  allows  multi-mission  capability  of  U.S.  ships  depending 
on  the  assigned  tasking.  British  Aerospace  Engineering  (BAE)  has  expanded  this  system 
by  adding  a  MK  25  Quad-Pack.  4 


This  canister  allows  a  cell  that  traditionally  holds  one  missile  to  store  four  ESSMs  in  the 
same  space.  This  modification  allows  a  ship  to  have  an  increased  self-defense  capability 
(British  Aerospace  Engineering  (BAE  Systems  2017). 


8  Cell  Module 


Figure  1.  MK  41  VLS  Module.  Each  eight-cell  module  is  housed  below  the 
ship’s  deck.  Each  cell  is  capable  of  holding  and  launching  a  range  of 
missile  types.  Source:  Batzler  (2011). 


Figure  2.  Arleigh-Burke  Class  Guided  Missile  Destroyers  (DDG  51-119).  These 
destroyers  have  up  to  96  cells  across  eight  installed  modules  in  the  aft 
part  of  the  ship  and  four  installed  modules  in  the  forward  part  of  the 
ship.  Source:  FAS  (2016b). 
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Figure  3.  Ticonderoga  Class  Guided  Missile  Cruisers  (CG  52-73).  These 
cruisers  have  up  to  122  cells  across  eight  installed  modules  in  the  aft 
part  of  the  ship  and  eight  installed  modules  in  the  forward  part  of  the 
ship.  Source:  FAS  (2000). 

5.  MK  57  Advanced  Vertical  Launching  System  (AVLS) 

The  introduction  of  the  Zumwalt  Guided  Missile  Destroyer  (DDG  1000),  shown 
in  Figure  4,  brought  with  it  a  new  vertical  launching  system  known  as  the  MK  57 
Advanced  Vertical  Launching  System  (AVLS)  (Raytheon  2012).  The  launcher  is 
designed  to  accommodate  all  existing  missile  types  while  allowing  room  for  extra  growth 
for  new  missile  types  and  sizes  (Raytheon  2012).  The  MK-57  launchers  are  contained 
and  protected  by  the  Peripheral  Vertical  Launching  System  (PVLS).  The  PVLS  modules 
distribute  the  ship’s  missile  launchers  in  separate  four-cell  launcher  compartments  along 
the  perimeter  of  the  ship’s  hull  (General  Dynamics  2013).  This  construction  makes  the 
launchers  and  contained  missiles  more  resistant  to  damage  incurred  from  enemy  attacks 
(Raytheon  2017b). 
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Figure  4.  Zumwalt  Class  Guided  Missile  Destroyers  (DDG  1000-1002).  These 
destroyers  have  up  to  80  cells  across  the  PVLS.  Image  Source: 
Raytheon  (2016  );  Caption  Source:  United  States  Navy  (2016a). 

B.  THESIS  CONTRIBUTION 

The  focus  of  this  research  is  to  develop  a  heuristic  solution  to  the  missile-to-VLS 
cell  optimization  problem  that  is  fast,  feasible,  and  nearly  optimal.  To  determine  how 
good  the  solutions  are,  we  can  compare  how  well  the  heuristic  model  compares  with 
optimal  results  from  Wiederholt’s  (2015)  VLP.  The  final  goal  of  this  study  is  to  socialize 
this  problem  and  solution  at  all  levels.  If  we  can  get  the  model  into  the  hands  of 
operational  planners  in  an  accessible  and  usable  form,  we  can  start  basing  our  loadouts  on 
actual  operational  plans  (OPLAN)  and  concept  plans  (CONPLAN).  This  will  give  fleet 
commanders  a  good  “what-if?”  tool  to  balance  the  potential  risk  to  mission 
accomplishment  with  the  cost  of  missile  movement. 

C.  ORGANIZATION 

Chapter  II  provides  a  literature  review  on  missile  loadout  and  mission  assignment 
problems  and  a  heuristic-based  solution  to  one  such  problem.  The  chapter  also  includes 
more  in-depth  coverage  of  Wiederholt’s  (2015)  work,  which  provides  the  starting  point 
for  this  thesis.  Chapter  III  provides  the  formulation  for  the  original  model.  Chapter  IV 
provides  the  methodology  behind  the  heuristics  of  the  Missile  Loadout  Planner  (MLP). 
Chapter  V  provides  a  comparison  of  the  optimal  solution  and  the  heuristic  solution  and 
provides  an  analysis  of  how  well  the  planner  functions.  Chapter  VI  provides 
recommendations  for  future  work. 
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II.  LITERATURE  REVIEW 


A.  INTRODUCTION 

Past  projects  for  Maritime  Operational  Planners  address  aspects  of  missile  loading 
and  mission  planning.  Staffs  could  use  these  decision  support  tools  in  their  planning 
efforts.  Previous  research  includes  a  simulation  model  to  fill  VLS  cells,  optimizing  ship 
to  mission  assignment  scheduling,  and  maritime  operational  planning. 

Next,  optimizing  TLAM  loadout  for  given  missions  went  through  all  of  the 
phases  desired  of  a  maritime  planning  tool.  It  started  as  an  optimization  model,  was 
reviewed  and  extended  as  an  optimization  model,  and  a  heuristic  was  tested.  The  model  is 
currently  being  used  to  decrease  the  workload  on  operational  planners. 

Finally,  we  examine  the  precursor  to  this  research,  namely,  Wiederholt’s  (2015) 
original  problem. 

B.  PREVIOUS  RESEARCH 

Jarek  (1994)  used  simulation  to  create  a  model  to  fill  VLS  cells  in  support  of 
surface-to-air  missile  (SAM)  requirements  as  required  for  Anti- Air  Warfare  (AAW) 
defense.  Two  cases  were  explored:  one  in  which  combat  air  patrol  (CAP)  was  able  to 
assist  against  the  attack,  and  one  without  such  assets. 

Dugan  (2007)  developed  the  Navy  Mission  Planner  (NMP).  This  integer 
programming  tool  can  assist  decision  makers  in  assigning  ships  to  missions  and 
identifying  dependencies  such  as  assigning  a  shotgun  ship  for  air  defense  along  with  a 
ship  assigned  to  Theater  Ballistic  Missile  Defense  (TBMD)  for  a  given  set  of  mission 
priorities.  This  model  is  run  on  a  fictional  case  set  on  the  Korean  Peninsula. 

Deleon  (2015)  used  an  integer  linear  program  to  develop  the  Navy  Operational 
Planner  (NOP)  to  help  decision  makers  with  maritime  operational  planning.  His  work 
explores  our  Navy’s  capability  to  accomplish  missions  as  quickly  as  possible.  While  his 
work  focuses  on  mine-hunting  countermeasures  (MCM)  scenarios,  the  analysis  can  be 
extended  to  other  missions.  The  model  considers  the  ability  of  a  ship  to  complete 
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multiple  missions  due  to  possible  dependencies  among  the  missions,  such  as  a  particular 
sequence  of  missions  that  must  be  adhered  to,  thus  saving  transition  time  and  achieving 
faster  mission  completions.  This  model  was  developed  to  advise  operational  planners  on 
more  efficient  mission  completion  than  the  traditional  manual  planning  method.  The 
traditional  method  quickly  becomes  time-consuming  and  prone  to  human  error  as  the 
number  of  missions  and  requirements  grow,  amplifying  the  need  for  a  computational  tool 
to  allocate  resources. 

C.  TLAM  OPTIMIZATION  SUPPORT  TOOL 

Kuykendall  (1998)  used  an  integer  programming  model  to  optimize  tomahawk 
land  attack  missiles  (TLAM)  to  strike  missions.  The  model  takes  into  consideration  the 
loadout  of  each  unit,  specific  tasking  and  geographic  location.  This  work  was  later 
generalized  by  Newman  et  al.  (201 1)  and  is  currently  in  use  by  the  U.S.  Navy. 

Newman  et  al.  (2011)  extended  conventional  weapon  assignment  problems  to 
develop  a  model  to  assist  operational  planners  in  the  Tomahawk  pre-designation  problem. 
This  model  advises  which  firing  platforms  should  be  selected,  and  on  each  platform  how 
many  and  which  type  of  Tomahawk  missiles  should  be  designated  to  meet  a  specific 
mission.  Since  firing  units  leave  port  with  a  preset,  fixed  loadout,  number,  and  variety  of 
Tomahawk  missiles,  a  tool  was  developed  to  best  meet  this  challenge.  This  is  a  mixed 
integer  linear  program  optimization  model  with  multiple  objective  functions.  One 
approach  presented  to  solve  the  multiple  objective  functions  is  to  scale  them  down.  A 
linear  weighted  sum  forms  a  linear  value  function  that  provides  an  accurate 
representation  of  the  problem  to  decision  makers.  Another  approach  is  hierarchical 
optimization  by  prioritizing  the  objective  functions  and  solving  the  objective  functions 
optimally  in  descending  priority  order.  The  solved  objectives  are  then  added  as 
constraints  requiring  some  fraction  of  their  achievement  by  the  following  problems. 
Alternatively,  a  fast  heuristic  algorithm  was  developed  to  find  good  solutions  quickly. 
Salmeron  (2002)  expands  the  details  of  this  heuristic.  This  model  was  then  tested  on  a 
series  of  problems  that  showed  the  heuristic  results  are  comparable  to  those  produced  by 
exact  mathematical  optimization.  All  the  computational  methods  proposed  by  these 
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researchers  performed  better  than  the  historic  method  of  a  given  Tomahawk  Strike 
Coordinator  (TSC)  manually  preparing  a  strike  plan. 

D.  ORIGINAL  VLS  LOADOUT  PROBLEM  AS  DESCRIBED  BY 

WIEDERHOLT 

Wiederholt  (2015)  developed  and  implemented  an  optimization  tool  known  as  the 
Vertical  Launch  System  Loadout  Planner  (VLP)  in  support  of  Seventh  Fleet  operational 
planners  in  the  Western  Pacific.  His  mixed  integer  linear  program  optimally  assigns 
missile  loadouts  to  ships  and  ships  to  missions  to  reduce  the  number  of  uncovered 
potential  missions.  This  model  is  implemented  in  the  General  Algebraic  Modeling 
System  (GAMS)  (2015).  His  work  considered  restrictions  such  as  a  limited  availability  of 
missile  type,  acceptable  substitutions  for  given  missiles,  priorities  on  missions,  missile 
compatibility  on  ships,  minimum  number  of  ships  required  for  mission  completion,  and 
changing  a  pre-existing  missile  loadout.  Violations  of  these  restrictions  lead  to  penalties 
in  the  model.  The  problem  was  formulated  by  looking  at  two  war  plans,  each  with  a  set  of 
missions  over  three  main  scenarios.  The  three  scenarios  covered  are  a  fixed  loadout, 
some  cells  fixed,  and  a  completely  unrestricted  scenario  where  the  VLP  is  allowed  to 
optimally  load  the  ships.  Wiederholt  (2015)  used  23  guided  missile  ships  with  the  multi¬ 
mission  module  of  the  MK41  Vertical  Launching  System  (VLS).  The  two  sets  of 
combatants  covered  are  Forward  Deployed  Naval  Forces  (FDNF)  and  West  Coast 
deploying  forces  on  two  six-month  deployment  cycles.  He  finds  that  as  the  fixed 
restrictions  decrease,  the  decision  tool  is  able  to  reduce  the  number  of  uncovered 
missions. 

E.  PYTHON  GENETIC  ALGORITHM  HEURISTIC 

The  basic  idea  and  sketch  for  the  genetic  algorithm  stemmed  from  a  heuristic 
formulated  by  CAPT  (ret.)  Jeff  Hyink  (2016)  and  executed  in  Python  using  Enthought 
Canopy  (2012).  This  simplified  version  of  the  GAMS  developed  solution  was  used  as  a 
proof  of  concept  of  the  viability  of  the  genetic  algorithm’s  ability  to  solve  the  problem. 
The  start  of  the  VBA  algorithm  rested  on  his  initial  heuristic  and  non-developed  areas 
were  instantiated  to  get  a  better  one-to-one  comparison  with  the  optimal  GAMS  solution. 
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III.  MODEL  FORMULATION  FOR  OPTIMIZING  VLS  MISSILE 

LOADOUT 


A.  INTRODUCTION 

Chapter  III  covers  the  formulation  for  the  Vertical  Launch  System  (VLS)  Loadout 
Planner  (VLP).  The  formulation  starts  by  building  the  indexes  starting  with  the  simplest 
then  accumulating  the  multi-dimensional  tuples.  Some  sets  are  created  in  the  course  of 
the  program  dynamically  to  only  represent  what  is  reflected  in  the  actual  data.  An  asterisk 
(*)  in  the  formulation  text  indicates  changes  in  wording  and  a  red  index  in  the 
formulation  indicates  a  change  of  index  from  the  original  model.  Constraint  (D3) 
proposes  a  modification  to  the  original  Wiederholt  formulation  to  ensure  that  if  a  mission 
is  not  committed  to,  then  there  is  no  ship  penalty  associated  with  that  shortage. 


B.  MODEL  FORMULATION  TO  OPTIMIZE  THE  MK  41  VERTICAL 
LAUNCH  SYSTEM:  VLP 


1.  Index  Use  [-Cardinality] 


W  G  W 
m  e  M 
d  g  D 
cgC 

S  E  S 
h  g  H 
y  eY 
teT 
reR 


warplan  [-10] 

missions  (alias  m)  [-10]  (e.g.,  TBMD  station) 
deployment  cycles  [-2] 

required  mission  ship  classes  (includes  class  “any”) 
[~6] 

individual  ships  [-25] 
home  ports  [-2] 

missile  types  (alias  y' ,  y  desired,  y'  committed)  [-8] 

type  of  mission  [-3] 

risk  level  (including  “high”)  [-2] 

class  of  ship  5 

type  of  mission  m 
risk  of  mission  m 


2.  Useful  Tuples 

(Those  marked  with  an  asterisk  are  derived  and  filtered  from  the  others 
defined  by  data.) 
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{w,m}  sWM  * 

missions  of  warp lan  w  [~  1 0] 

{w,d,m}  g  WDM 

warplan-mission-cycle  triples  [10x10x2] 

{m,  c}  g  MC 

mission  m  can  be  completed  by  ship  class  c 

{w,d, m,  5}  g  WDMS  * 

warplan-mission-cycle-ship  4-tuples  [10x10x2x25] 

KjJgSF 

ship  5  cannot  accommodate  missile  type  y 

{s,d}eSD 

ship  5  deployment  cycles 

{w,  d,  m,  v}  g  WDMY  * 

warplan-mission-cycle-missile  4-tuples  [10x10x2x10] 

{m,m'}GMM'* 

missions  m  and  m '  are  mutually  exclusive 

(e-g 

{m,y,y'}  eMYY'  * 

missile  type  y  can  be  substituted  for  type  y’ 

{w, d,  m,  s,  y}  g  WDMSY  * 

5-tuple  for  missile  requirements,  or  loading 

{ w,  d,  m,  s,  y,  y ' }  g  WDMSYY '  * 

6-tuple  for  missile  loading  with  substitutions 

3.  Given  Data  [Units] 

priority  m 

priority  of  mission  m  [penalty] 

ships _reqm 

ships  required  by  mission  m  [ships] 

ships_short _penm 

ship  shortfall  penalty  for  mission  m  [penalty /ship] 

missiles  _desiredmy 

desired  type  y  missiles  *on  each  ship  for  mission  m 

[missiles] 

missiles_minimumm 

minimum  missiles  *on  each  ship  of  type  y  for  mission 

m  [missiles] 

missile_short _penm  y 

missile  shortfall  penalty  for  mission  m,  type  y  [penalty/ 
missile] 

vls_cellss 

number  of  VLS  cells  on  ship  s  [cells] 

missile  inventory 

number  of  type  y  missiles  in  inventory  [missiles] 

missiles  _per_cellx 

number  of  type  y  missiles  in  a  VLS  cell  [missiles  per 
cell] 

risk_missileJoads 

ship  s,  type  y  missiles  carried  in  addition  to  mission 
load  [missiles] 

under _peny,over _penv 

penalty  for  disproportionate  spread  of  missile  type  y 
among  ships  carrying  these  for  each  mission  [penalty] 

min_missile_load  s 

ship  5  type  y  missiles  carried  in  addition  to  mission  load 

[missiles] 

alt_missile  _penm  , 

penalty  for  substituting  type  y' for  y  in  mission  m 
[penalty/missile] 
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loadouts  y,  ship  s  load  of  missile  type  y'  prior  to  optimization 

[missiles] 

change _pen  penalty  for  adjusting  prior  loadout  [penalty/missile] 


4.  Decision  Variables  [Units] 


ASSIGN 


w,d,m,s 


MISSION. 


COMMIT. 


w,d  ,m,s,y,y' 


MISSILE  SLACK 


w,d,m,y 


SHIPS  SHORT 


w,d  ,m 


MISSILES  SHORT. 


w,d,m,y' 


LOAD 


assign  ship  s  to  warplan  w  deployment  cycle  d  mission 
m  [binary] 

plan  w,  cycle  d,  mission  m  commitment  [binary] 
plan  w,  cycle  d,  mission  m,  want  type  y,  commit  type 
/  [missiles] 

plan  w,  cycle  d,  mission  m,  type  y  missiles  short  of 
desired  number  [missiles] 

plan  w,  cycle  d,  mission  m,  elastic  variable  for  ship 
shortages  on  mission  [ships] 

plan  w,  cycle  d,  mission  m,  elastic  variable  for  type  y 
missile  shortages  on  missions  [missiles] 
ship  5  load  of  missile  type  y’  [missiles] 


RISK_MISSILESsy, 
UNDERwdmsy ,  O  VERt 
CHANGE^, 


DEPLOY s 
DEPLOY _WARws 


carried  by  ship  assigned  high-risk  mission(s)  [missiles] 
elastic  variable  for  inequitable  missile  loads  [fraction] 
number  of  y’  missiles  changed  in  VLS  cells  of  ship  s 
[missiles] 

indicator  that  ship  5  is  deployed  [binary] 

indicator  that  ship  5  is  deployed  in  war  plan  w  [binary] 
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5. 


Formulation 


min  - 

ASSIGN, 
MISSION, 
COMMIT, 
MISSILE_SLACK, 
SHIPS_SHORT, 
MISSILES _SHORT, 
LOAD, 

RISK_MISSILES, 

UNDER.OVER 


X  priority  MISSION w4  t 

{ w,d,m  }eWDM 


s.t. 


+  ^  altjnissile _penm  ,  COMMIT wdms  . 

{ w,d,m,s,y,y'}eWDMSYY ' 

+  ^  ships_short _penSHIPS_SHORTwJ 

{ w,d ,m  }eWDM 

+  ^  missiles jshort _penm  MISSILES_SHORT  d  , 

{ w,d,m,y'}eWDMSY 

+  ^  under _pen  UNDER  d  ^  ^ 

{ w,d ,m,s,y  }eWDMSY 

+  ^  over  jjen  O  VER  d  m 

{ w,d  ,m,s,y  }eWDMSY 

+  ^  change _penCHANGE  . 

[s.y'HSY 


ASSIGN  ,  +  ASSIGN  ,  <1 


MISSION  ,  >  ASSIGN 


V{w,d,m,s }  e  WDMS, 
{w,d,m',s}  eWDMS  I 
{m,m'}  e  MM' 
\/{w,d,m,s}  e  WDMS 


X  ASSIGN  w 

s\{w,d,m,s}<=WMDS 


+ SHIPS _SHORT  d  m  =  ships_req 


MISSION w  d  m  m]  e  WDM 


RISK  MISSILES 

s,y 


>  (min  missile  load  , 

■  —  —  s,y 

+  riskjnissilejoad  s  y,  I,  =,  hish.)  ASSIGN  nd  \/{w,d,m,s,y'}  gWDMSY 

Y  COMMIT  „ 

w,d  ,m,s  ,y,y 

y'\{w,d,m,s,y,y'}<=WDMSYY' 


<  missiles_desired  ASSIGN w dmi  \/{w,d,  m,  s,  y }  e  WDMSY 

LOAD  ,  >  V  COMMIT  , 

s,y  w,d  ,m,s  ,y,y 

{d ,m,y}\{  w,d ,m,s ,y  ,y' }gWDMSYY' 

+  RISK  MISSILES  ,  Vw  e  W,{,v,  v'}  <£  SY 


(DO) 


(Dl) 
( D2) 


(D3) 


(D4) 


( D5) 


(D6) 
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COMMIT  d  m 

{ s,y  '}l{  w,d,m,s,y,y  '}eWDMSYY  ' 


+MISSILES_SHORTwdmy  +  MISSILE_SLACKwd 


>  missiles  desired  MISSION  , 

—  m,y  w,d  ,m 

V{w,d,m,y}  g  WDA/7 

(D7) 

x-  1 

/  - LOAD  <  vls_cells 

[s,r]tsr  missiles _per_cellx, 

Vi  g  5 

(D8) 

^  LOAD  v,  <  missile -inventory  , 

{s,y')iSY 

v_v 1  g  y 

(D9) 

V  COMMIT  +  UNDER 

w,d  ,m,s,y,y  w,d,m,s,y 

y  'l{  w,d  ,tn  ,s ,  y ,  y  ’ )  eWDMSyy  ' 

OVER  d  m  s  v 

=  (missiles_desiredm  /  ships_reqm) ASSIGN w d  m  s 

V{w, <7, /w,  i,  y }  g  WDMSY 

(DIO) 

CHANGE  v  >  +( LOAD  -  loadout ,  ,) 

V{s,y'}tSY\  loadout  >0 

"  y 

yl{5,y}gSF 

(Dll) 

CHANGE  >  (LOAD  loadout  I 

v{*,y'}esn  s  ,oaiw  *> 

yl{5,y}gSr 

(D12) 

DEPLOY ;  >  ASSIGN wdmi 

V{w,</,m,i}  gWDMS 

(D13) 

DEPLOY  <  J]  ASSIGN  „  d  m  , 

/nl{vi',rf,m,s}elVDMS 

V{w,d,s]  <=WDS 

(D14) 

DEPLOY_WAR  >  ASSIGN wd  m 

V[w,d,m,s]  eWDMS 

(D15) 

DEPLO  Y_WAR  _  <  ^  ASSIGN  ,  m  i 

ml{  w,d,/n,s)eWDM5 

V{w,d,s}  eWDS 

(D16) 

ASSIGN wd  ms  g  {0,1} 

\/{w,d,m,s }  g  WDMS 

MISSION  ,  g  {0,1} 

V{w,d,m]  g  WDM 

COMMIT  ,  g  Z+ 

\/{w,d,m,s,  y,  y '} 

g  WDMSYYP 

0  <  MISSILE  SLACK  , 

—  w,d,m,y 

<  missiles  desired  -  missiles  minimum 

—  m,y  —  l 

n  v  V{  w,  d,m,y}  g  WDMY 

SHIPS_SHORT  >  0 

\/{w,d,m }  g  WDM 

MISSILES  SHORT  >  0 

w,d,m,y 

V{w,d,m,  y}  g  WDMY 

LOAD  . >  0 

V{i,  y}  g  57 

RISK  M ISSILES  >  0 

V{i,  y}  g  57 

UNDER  d  m  sy„ OVER  dmsy,>0 

V{w,  i,  y '} 

g  WDMSY 

CHANGE  > 0 

*>y 

{i,  y '}  g  57 

DEPLOY  >  0 

Vi  g  5 

DEPLOY _WAR  s  >  0 

Vw  G  W,  i  G  5 

(D17) 
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6. 


Discussion 


This  optimization  model  provides  the  best  single  VLS  loadout  for  each  ship.  Each 
ship  is  given  one  loadout  regardless  whether  she  is  a  deployer  making  one  deployment 
cycle  or  a  FDNF  ship  making  two  cycles.  The  solution  also  advises  the  best  ship-to- 
mission  pairing  for  some  number  of  warplans.  The  loadouts  provide  the  best  solution  to 
be  prepared  regardless  of  which  warplan  is  required.  Planners  retain  control  to  manually 
set  a  ship’s  loadout.  The  optimization  model  will  account  for  this  fixed  loadout  and 
assign  the  remaining  ship  loadouts  and  mission  assignments  to  best  cover  the  required 
missions.  Wiederholt  describes  the  objective  function  and  constraints  of  this  model: 

1.  The  objective  (DO)  accounts  rewards  for  prioritized  mission 
accomplishment  and  deducts  penalties  for  violating  policies  that  cannot  be 
satisfied.  A  number  of  these  penalties  result  from  optional  model  features. 

2.  Each  constraint  (Dl)  restricts  a  ship  from  performing  mutually-exclusive 
missions. 

3.  Each  constraint  (D2)  signals  a  mission  accomplishment  if  any  ship  is 
assigned  to  this  mission. 

4.  Each  constraint  (D3)  provides  the  required  number  of  ships  for  a  mission, 
or  accounts  for  any  shortfall. 

5.  Each  constraint  (D4)  reckons  whether  a  ship  needs  extra  defensive 
missiles  due  to  the  risk  level  of  missions  assigned  to  it.  This  is  later 
referred  to  as  a  “risk  load”  or  “defensive  load.” 

6.  Each  constraint  (D5)  commits  a  number  of  a  required  missile  type,  or  an 
acceptable  substitute  type  to  fulfill  an  assigned  mission. 

7.  Each  constraint  (D6)  reckons  the  number  of  missiles  of  some  type  that  are 
to  be  loaded  on  a  ship. 

8.  Each  constraint  (D7)  reckons  whether  the  required  number  of  missiles  has 
been  loaded,  or  accounts  for  a  shortfall. 

9.  Each  constraint  (D8)  limits  the  number  of  missiles  that  can  be  loaded  into 
the  vertical  launching  system  of  a  ship. 

10.  Each  constraint  (D9)  limits  the  number  of  a  type  of  missile  to  the  total  in 
inventory. 

11.  Each  (optional)  constraint  (DIO)  requires  that  a  type  of  missile  be  loaded 
proportionately  on  each  ship  participating  in  a  mission. 
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12.  Each  constraint  (Dll)  and  its  pair  (D12)  (optionally)  reckon  the  positive 
difference  between  a  pre-existing  VLS  loadout  and  the  one  being 
prescribed  by  the  model.  This  positive  difference  is  penalized  in  the 
objective  function  in  order  to  reduce  unnecessary  “turbulence”  between 
legacy  loadouts  and  their  optimal  revisions,  but  could  just  as  well  be 
limited  numerically  by  ship  and  by  missile  type  if  it  is  anticipated  that 
there  will  be  limited  pier-side  time  to  make  changes. 

13.  Constraints  (D13-14)  are,  together,  optional.  Each  constraint  (D13)  sets  an 
indicator  that  a  ship  has  been  assigned  a  mission  in  some  deployment 
cycle  of  some  war  plan.  Each  constraint  (D 1 4)  assures  that  a  deployed 
ship  is  assigned  at  least  one  mission  in  each  deployment  cycle  of  each  war 
plan. 

14.  Constraints  (D15-16)  are,  together,  optional,  and  are  subsumed  if 
constraints  (D13-14)  are  invoked.  Each  constraint  (D 15)  sets  an  indicator 
that  a  ship  has  been  assigned  a  mission  in  some  deployment  cycle  of  some 
war  plan.  Each  constraint  (D16)  assures  that  a  deployed  ship  is  assigned  at 
least  one  mission  in  each  deployment  cycle  of  each  war  plan  to  which  it 
has  been  assigned  a  mission. 

15.  Constraint  (D17)  defines  decision  variable  domains.  (Wiederholt  2015) 
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IV.  DATA  AND  METHODOLOGY 


A.  METHODS  FOR  HEURISTIC  ALGORITHMS 

Goldberg  (1989,  p.  6)  emphasizes  that  optimization  is  as  much  about  the  process 
of  improvement  as  it  is  finding  the  optimal  point.  The  appeal  of  an  exact  optimal  solution 
is  the  amount  of  certainty  that  can  be  provided  under  given  model  assumptions  and 
objective  function  (Rardin  1998,  p.  16).  Reeves  (1993,  p.  6)  defines  a  heuristic  as  a 
method  that  seeks  near-optimal  solutions  at  a  reasonable  computational  cost.  Effective 
heuristic  algorithms  enable  a  technique  to  produce  fast  approximate  solutions  that  are  not 
guaranteed  to  be  optimal  but  are  feasible  to  the  problem. 

Rardin  describes  methods  for  heuristic  algorithms  to  find  good  feasible  solutions 
that  are  approximately  optimal  (1998  pp.  15-16).  For  consistency,  this  thesis  will  use 
Rardin’s  notation  to  the  greatest  extent  possible.  A  solution  will  be  defined  as  a  choice 
for  the  decision  variables  denoted  by  the  vector,  x.  A  component  is  a  scalar  member  of 
that  vector,  x(i),  where  i  denotes  the  index  in  the  solution  vector  of  the  scalar.  The  first 
solution  visited  by  the  heuristic  is  labeled  x(<)),  followed  by  x(1),  and  continues  until  a 
user-defined  limit  of  iterations  (Rardin  1998,  pp.  78-79).  Generally,  improving  searches 
over  discrete  variables  (integer  or  binary)  are  defined  in  a  specified  neighborhood  around 
the  current  solution,  x(t),  for  each  iteration  t.  The  process  continues  until  a  user-defined 
iteration  limit  has  elapsed.  This  neighborhood  contains  a  set  of  allowable  nearby 
solutions  that  can  be  reached  from  the  current  solution  by  a  simple  operation  (Reeves 
1993,  p.  5).  This  set,  known  as  a  move  set ,  contains  the  current  solution  and  all  solutions 
within  a  small  adjustment  of  discrete  values  in  the  area  around  it.  Each  iteration  of  the 
algorithm  checks  for  feasibility  and  a  superior  objective  value  within  this  discrete 
neighborhood.  Once  no  such  improvement  exists,  a  local  optimum  is  attained.  Figure  5 
depicts  a  local  versus  global  optimum. 
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Figure  5.  Example  Three-Dimensional  Surface  with  Two  Local  Maxima.  Only 
one  of  these  maxima  is  global.  A  constructive  heuristic  risks 
converging  to  the  lower  optima. 

B.  GENETIC  ALGORITHMS 

Genetic  search  algorithms  are  based  on  the  biological  principles  of  natural 
selection.  This  class  of  algorithms,  known  as  evolutionary  algorithms,  draws  an  analogy 
between  genetic  structures  and  complex  vector  components  of  an  optimization  problem 
(Reeves  1993,  p.  152).  Holland  (1975)  introduced  the  notion  of  a  progressive 
modification  of  genetic  structures  using  a  set  of  operators.  These  structures  adapt  using 
these  operators  and  are  evaluated  based  on  how  well,  or  fit  (healthy),  they  are  based  on 
some  measure  of  performance.  The  method  of  a  genetic  algorithm  blends  a  survival  of 
the  fittest  mentality  and  a  human-influenced  structured  random  search.  This  blend,  when 
properly  balanced,  provides  a  robust  and  efficient  method  for  finding  improving  solutions 
(Goldberg  1989,  p.  2).  Rardin  (1998,  p.  696)  describes  the  elitist  strategy  of  genetic 
algorithms.  This  strategy  forms  each  new  generation  as  a  mix  of  the  elite,  or  best, 
solutions  carried  over  from  previous  solutions,  arbitrary  solutions,  known  as  immigrants, 
to  increase  diversity  in  the  population,  and  children  of  crossover  operations  in  the 
previous  population.  This  diversity  enables  the  algorithm  to  search  new  parts  of  the 
feasible  region  in  a  search  for  new  solutions.  Some  of  these  solutions  may  even  have  poor 
objective  values. 

The  first  step  in  this  algorithm  is  to  define  the  genetic  representation  of  the 
problem.  This  representation  is  just  a  numeric  vector  of  length  n  known  as  the 
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chromosome  that  represents  possible  solutions,  x(t),  to  an  optimization  problem.  Positions 
in  the  bitstring  are  the  loci  of  the  chromosomes.  The  variable  at  a  given  locus  is  a  gene, 
and  a  gene’s  value  an  allele.  The  set  of  chromosomes  is  a  genotype.  This  genotype 
defines  an  individual  .phenotype.,  with  certain  fitness  (Aarts  and  Lenstra  1997,  p.  141). 

The  next  step  is  to  create  a  random  initial  population  of  possible  chromosomes,  or 
solutions,  denoted  P(0).  Chromosomes  are  evaluated  according  to  a  fitness,  or  objective, 
function.  This  fitness  function  is  some  measure  of  goodness.  The  most  fit  gene  pool  will 
survive  to  the  next  generation’s  population  and  continue  to  produce  better  solutions 
(Ragsdale  2008,  p.  381).  Figure  6  gives  an  example  of  an  initial  population  and  genetic 
representation 

Consider  the  following  simple  optimization  as  an  example: 

max  xl  +  x2  +  x3  +  x4  +  xs 

a: 

s.t.  *,.{0,1}  Vi 


INITIAL  POPULATION 


Figure  6.  Example  Genetic  Representation  of  an  Optimization  Model.  This 
example  shows  a  set  of  randomly  generated  gene  values.  Each 
chromosome  and  its  respective  evaluated  fitness  value  are  listed  in  a 
random  generation  of  the  initial  population. 
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Goldberg  (1989,  p.  10)  gives  three  simple  operations  for  forming  the  next 
generation:  reproduction ,  crossover,  and  mutation.  Reproduction  is  simply  copying  the 
top  performing  individual  strings  by  objective  value  into  the  next  generation.  Each 
generation  is  known  as  the  population  at  iteration  t,  denoted  P(t).  After  reproduction, 
crossover  operations  occur  by  selecting  N  ( N  even)  chromosome  vectors  from  P(t)  as 
parents.  Parents  are  paired  to  form  N/2  pairs.  Each  pair  of  parents  is  crossed  by  selecting 
a  random  integer  position  along  the  chromosome  string  and  splitting  both  parent  vectors 
at  that  position.  Children  are  formed  by  recombining  the  part  above  the  position  of  one 
parent  solution  with  the  part  below  the  position  of  the  other  parent.  The  same  is  done  for 
the  remaining  parent.  Both  children  become  members  of  the  new  population.  It  is 
important  to  check  solutions  after  crossover  to  ensure  the  resulting  children  remain 
feasible.  This  does  not,  however,  guarantee  improvement  in  the  objective  function 
(Rardin  1998,  p.  695).  The  mutation  operation  is  the  random  replacement  of  values  in  a 
solution  vector.  Figure  7  shows  an  example  crossover  operation  and  Figure  8  shows  an 
example  mutation  operation.  Figure  9  shows  the  resulting  fitness  values  after  all 
operations  are  complete  and  Figure  10  shows  the  new  population  formed  by  keeping  the 
top  performers. 
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Figure  7.  Crossover  Operations  on  Parents  PI  and  P2.  Chromosomes  PI  and  P2 
are  chosen  as  parent  vectors.  Crossover  operations  select  a  random 
position  in  the  parent  vectors  and  exchange  the  latter  parts  to  form  two 
children  chromosomes  Cl  and  C2  that  are  then  added  to  the  next 
generation  of  solution  vectors. 
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MUTATION 
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Figure  8.  Mutation  Operations  on  a  Chromosome.  Mutation  operations  choose 
one  or  two  random  alleles  to  flip  in  a  chromosome  vector.  The 
resulting  vector  is  then  introduced  into  the  set  of  most  fit  candidates 
for  the  next  generation. 


CROSSOVER  AND  MUTATION 


MUTATION 


CROSSOVER 


Figure  9.  Subsequent  Generations  are  Based  on  Fitness  Values.  Fitness  values 
are  calculated  for  the  child  chromosomes  formed  from  crossover 
operations  and  allele  flips  from  the  mutation  operations.  The  top 
performers  are  carried  on  to  the  next  generation. 
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NEW  POPULATION 
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Figure  10.  Top  Performers  Form  the  New  Population.  The  resulting  population 
after  one  iteration  of  a  genetic  algorithm  is  formed  by  taking  the  top 
performers  in  terms  of  fitness  in  the  previous  generation  from 
reproduction,  children  of  crossover  operations  and  allele  flips  from 
mutation  operations. 

This  process  of  crossover  and  mutation  to  form  new  generations  of  populations 
continues  until  a  user-defined  iteration  limit  has  elapsed.  The  best  solution  at  that  time  is 
taken  as  the  optimal  solution.  There  is  no  restricted  neighborhood,  so  there  is  little  need 
to  worry  about  getting  stuck  in  local  optima. 

1.  Abstract  Algorithm 

Separate  comma  separated  value  (CSV)  files  provide  all  of  the  input  data  sets 
required  for  the  model.  The  data  is  read  and  compiled  into  one  Microsoft  Excel  (2010) 
workbook.  Next,  all  of  the  data  and  index  sets  are  stored  internal  to  VBA  (2010)  as 
collections.  The  model  sets  are  built  starting  from  the  simplest  one-dimensional  sets.  An 
example  of  one  such  set  from  the  model  formulation  in  Chapter  III  is  the  set  of  warplans, 
w.  Next,  the  given  fixed  multi-dimensional  sets  are  built  and  stored.  The  set  of  missions 
m  that  can  be  completed  by  ship  class  c,  denoted  {m,c}  in  the  model  formulation  from 
Chapter  III  is  an  example  of  one  fixed  multi-dimensional  set.  Finally,  the  series  of 
dynamic  multi-dimensional  sets  based  on  tuples  of  fixed  sets  are  built  and  stored.  An 
example  of  dynamic  tuple  from  Chapter  III  is  the  set  of  missions  m  of  warplan  w,  denoted 
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{w,m}.  This  set  is  defined  by  the  data;  for  example,  each  warplan  may  or  may  not  have 
every  mission  as  a  requirement.  Constructing  the  sets  in  this  fashion  keeps  set  dimensions 
as  small  as  possible  which  leads  to  faster  computation  times  because  we  only  need  to 
explore  those  solutions  that  are  actually  feasible  and  not  every  possible  Cartesian  product 
of  the  component  sets.  The  data  is  stored  in  VBA  dictionaries  to  access  values  for 
penalties  and  feasibility  checks.  VBA  user-defined  types  are  created  to  identify  the  base 
set  of  missions  and  base  set  of  ships  and  their  associated  attributes  such  as  number  of 
ships  and  missiles  required  for  a  mission  and  ship  class  and  current  loadout  for  a  ship. 
This  base  remains  unchanged  while  the  genetic  algorithm  uses  simple  numeric  array 
values  to  reference  this  information. 

Each  chromosome  is  defined  by  a  loadplan.  The  loadplan  describes  a  full 
mapping  of  missiles  to  ships  and  ships  to  missions.  Each  loadplan  is  filled  with  decision 
variables,  called  genes,  defined  as  an  array  of  random  ship  assignments.  This  array  is  the 
same  size  as  the  total  number  of  ships  required  for  all  missions.  During  this  ship 
allocation,  feasibility  checks  are  also  performed.  Infeasible  solutions  will  be  rejected 
from  the  initial  population. 

Each  feasible  ship  assignment  is  then  mapped  to  a  set  of  missions  until  that 
mission’s  ship  requirement  is  filled.  Each  ship  assigned  to  missions  is  then  topped  off 
from  the  pier  inventory  to  best  meet  mission  minimum  missile  requirements  as  long  as 
there  are  empty  VLS  cells  available.  The  top-off  sequence  is  a  randomly  chosen  order  of 
the  ships  assigned.  Each  ship,  in  turn,  greedily  fills  all  of  her  empty  cells  as  long  as  there 
is  a  missile  shortage  for  an  associated  mission. 

After  the  ship  assignment  and  missile  top-off  stage,  a  fitness  array  for  each 
chromosome  is  calculated  using  the  penalties  from  the  same  Wiederholt  objective 
function. 

The  Genetic  Algorithm  Heuristic  (Rardin,  1998,  pp.  694-696): 

1.  The  first  generation  is  a  random  generation  of  N  constructed  solutions,  or 
loadplan  schemes.  Each  initial  chromosome  is  formed  by  iterating  over 
each  mission,  and  for  each  ship  requirement  in  the  mission,  randomly 
assign  a  compatible  ship.  The  ship  is  rejected  if  it  is  already  assigned  to 
the  mission  or  if  the  ship  is  already  assigned  to  another  mission 
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incompatible  with  the  mission.  This  process  or  random  ship  assignment  is 
repeated  until  the  mission  has  a  sufficient  number  of  ships. 


2.  After  the  ship  assignment  step,  use  a  greedy  loading  heuristic.  It  begins 
with  a  fresh  pier  inventory  to  top  off  the  ships  as  if  each  ship  were  empty. 
The  sequence  of  top  off  is  generated  by  a  random  order  of  the  assigned 
ships  and  a  random  order  of  missile  loading  on  each  ship  and  filled  in  a 
greedy  fashion.  Iterate  through  the  random  sequence  of  ships.  Proceed  to 
load  as  long  as  there  are  available  cells  on  each  ship  and  remaining  pier 
inventory.  For  each  defensive  missile  requirement  based  on  the  mission 
assigned  fill  the  defensive  load  onto  each  ship.  Once  each  ship  has  its 
defensive  load,  for  each  ship  assigned  to  a  mission,  try  to  fill  each  ship  to 
the  minimum  missile  requirement  for  each  missile  and  each  mission.  After 
defensive  and  mission  required  missiles  are  loaded,  attempt  to  top  off  each 
ship  to  its  original  loadout  as  long  as  cells  and  inventory  remains.  Capture 
the  new  change  from  the  old  loadout  to  the  new  loadout  for  the  fitness 
evaluation. 

3.  Calculate  the  fitness  value  of  each  of  these  loadplans  using  the  objective 
function  from  the  original  Wiederholt  (2015)  formulation.  All  solutions  in 
the  initial  population  are  generated  such  that  they  are  all  feasible  to  the 
problem.  Missions  that  are  not  assigned  the  minimum  number  of  missiles 
are  marked  as  incomplete. 

4.  Perform  a  sequence  of  mutation  procedures.  A  simple  mutation  will 
randomly  select  a  loadplan  from  the  initial  population  to  mutate.  From  this 
loadplan,  randomly  select  one  ship.  Perform  a  check  to  see  if  the  ship 
assignment  is  fixed.  If  it  is  fixed,  select  another  ship.  Once  a  ship  valid  to 
conduct  a  swap  is  selected,  perform  checks  on  each  mission  to  see  if  the 
selected  ship  can  fill  that  mission.  Also,  perform  checks  to  see  which  other 
ship  in  the  remaining  chromosome  is  eligible  to  fill  selected  ship’s 
mission.  Perform  the  loading  heuristic  on  each  eligible  swap,  calculate  the 
fitness,  and  keep  the  solution  if  the  prospective  loadplan  fitness  value  is 
less  than  the  old  loadplan  fitness  value. 

5.  The  next  mutation  performed  is  similar  to  the  first.  An  eligible  random 
ship  is  selected.  The  selected  ship  is  then  removed  from  the  chromosome 
and  replaced  by  an  eligible  ship  outside  the  chromosome  from  the  set  of 
ships.  Every  possible  replacement  for  the  selected  ship  is  tried  followed  by 
the  loading  heuristic  and  fitness  calculation.  The  loadplan  with  the  least 
fitness  is  kept  after  each  trial. 

6.  After  some  number  of  mutations,  perform  natural  selection.  This  evaluates 
the  fitness  of  two  randomly  selected  loadplans.  The  weaker,  or  less  fit, 
loadplan  is  replaced  with  the  more  superior.  This  results  in  two  copies  of 
the  same  more-fit  loadplan  that  will  then  be  part  of  the  next  mutation 
sequence. 
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7.  Continue  to  perform  mutations  until  some  user-defined  limit. 

8.  At  the  end  of  the  genetic  algorithm,  the  best  loadplan  is  pulled  out  and  is 
run  through  the  loading  heuristic  many  times  to  try  to  improve  the  ship 
loadouts,  and  thus,  refine  the  fitness. 
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V.  ANALYSIS  AND  RESULTS 


A.  INTRODUCTION 

Widerholt  (2015)  scenarios  serve  as  a  benchmark  for  the  Missile  Loadout  Planner 
genetic  algorithm  heuristic.  Common  warplan  data  and  penalties  will  be  used  to  evaluate 
three  scenarios.  The  first  scenario  is  tested  with  a  fixed  loadout  for  each  ship  and 
measured  according  to  its  ability  to  fulfill  mission  requirements  set  forth  by  Commander 
Seventh  Fleet  (C7F).  The  second  scenario  will  allow  the  Missile  Loadout  Planner  to 
choose  missile  loadouts  for  West  Coast  ships  and  keep  the  Forward  Deployed  Naval 

Forces  (FDNF)  ship  loadouts  fixed.  The  final  scenario  will  allow  the  Missile  Loadout 

Planner  to  choose  a  missile  loadout  for  each  ship  and  make  mission  assignments. 

B.  COMPUTATION  DETAILS 

An  initial  population  for  each  scenario  of  200  complete  loadplans  will  be 
generated.  Following  this  initial  generation,  a  repeated  sequence  of  mutations  followed 
by  natural  selection  leads  to  identification  of  the  most  fit  gene  evaluated  according  to  the 
objective  function  as  described  by  the  Wiederholt  (2015)  formulation  discussed  in 
Chapter  III.  Once  the  data  is  input  by  the  operational  planner,  the  tool  takes  an  initial 
overhead  time  of  roughly  20  minutes  to  internally  build  all  of  the  sets  and  data  required 
for  the  algorithm. 

C.  WARPLAN  SCENARIO  DATA 

The  data  is  a  replication  of  the  model  used  by  Wiederholt.  The  scenario  consists 
of  two  warplans  consisting  of  22  missions  and  30  missions  respectively  over  a  period  of 
two  deployment  cycles.  There  are  23  warships  available  for  use  that  may  or  may  not  be 
available  in  a  given  deployment  cycle  and  may  or  may  not  be  eligible  for  a  given 
mission.  Table  1  outlines  ship,  mission,  and  missile  abbreviations  that  will  be  used  going 
forward.  Table  2  denotes  warplan-mission  priorities,  or  the  score  for  committing  to  a 
mission.  Table  3  denotes  war-plan  mission  risk  types,  the  minimum  ships  required  to 
fulfill  the  mission,  and  a  penalty  for  each  ship  short  of  the  requirement.  Table  4  provides 
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the  reference  for  the  overall  missile  inventory  to  be  distributed  across  the  ships  and 
excess  to  remain  on  the  pier.  Table  5  shows  missile  requirements  for  a  given  mission.  A 
mission  desires  some  number  of  missiles  but  can  be  satisfied  with  a  minimum  number  of 
missiles.  Missile  loads  less  than  desired  result  in  some  penalty  for  assumed  warfighting 
risk  for  each  missile  shortfall.  Table  6  outlines  potential  missile  substitution  for  a  desired 
mission-missile  requirement.  This  less  desirable  missile  will  incur  some  penalty  on  the 
fitness  evaluation  for  a  loadplan.  The  missile  shortage  penalty  in  Table  5  is  high  relative 
to  all  other  penalties.  It  is  vital  that  if  we  commit  to  a  mission,  we  have  the  missiles 
onboard  to  accomplish  that  mission. 


Table  1.  Missile  Loadout  Planner  Problem  Input:  Abbreviations  and  Acronyms. 

These  abbreviations  and  acronyms  are  used  in  the  original  problem 
and  carry  through  the  genetic  algorithm  model.  Source:  Wiederholt 
(2015). 


US  Warships 

Designation 

Ticonderoga  Guided  Missile  Cruiser 

CG  (52-73) 

Arleigh-Burke  Guided  Missile  Destroyer 

DDG  (51-106) 

Zumwalt  Guided  Missile  Destroyer 

DDG  1000 

Missions 

Abbreviation 

Theater  Ballistic  Missile  Defense 

TBMD 

Escort 

Escort 

Surface  Action  Group 

SAG 

Strike 

STRIKE 

Missiles 

Designation 

Associated  Mission(s) 

Tomahawk  (1-3) 

TLAM 

STRIKE 

Standard  Missile  2  Medium  Range 

SM2  MR 

Escort/SAG 

Standard  Missile  2  Extended  Range 

SM2  ER 

Escort/SAG/TBMD 

Standard  Missile  3 

SM3 

TBMD 

Standard  Missile  6 

SM6 

ESCORT/SAG/TBMD 

Anti-Submarine  Rocket 

ASROC 

ESCORT/SAG 
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Table  2.  Description  of  Warplan-Mission  Priorities.  A  and  B  represent  two- 
fictitious  warplans.  Warplan  A  has  11  missions  and  Warplan  B  has  15 
missions.  Source:  Wiederholt  (2015). 


Warplan 

Mission 

Priority 

A 

TBMD 

80 

A 

STRIKE  1 

40 

A 

STRIKE  2 

40 

A 

STRIKE  3 

40 

A 

Escort  1 

30 

A 

Escort  2 

30 

A 

Escort  3 

30 

A 

Escort  4 

30 

A 

SAG  1 

80 

A 

SAG  2 

60 

A 

SAG  3 

60 

B 

TBMD  1 

70 

B 

TBMD  2 

70 

B 

TBMD  3 

60 

B 

Escort  1 

50 

B 

Escort  2 

50 

B 

Escort  3 

50 

B 

Escort  4 

50 

B 

Escort  5 

50 

B 

Escort  6 

50 

B 

SAG  1 

80 

B 

SAG  2 

50 

B 

SAG  3 

50 

B 

SAG  4 

80 

B 

SAG  5 

50 

B 

SAG  6 

50 
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Table  3.  Description  of  Warplan-Mission  Attributes.  These  attributes  include 
risk  type,  minimum  ship  requirement,  and  ship  shortage  penalty. 
Source:  Wiederholt  (2015). 


Warplan 

Mission 

Risk  Types 

Minimum 

Ships 

Ship  Shortage  Penalty 

A 

TBMD 

Very  High 

2 

100 

A 

STRIKE  1 

Medium 

1 

100 

A 

STRIKE  2 

Medium 

1 

50 

A 

STRIKE  3 

Medium 

1 

50 

A 

Escort  1 

Medium  Low 

1 

100 

A 

Escort  2 

Medium  Low 

1 

100 

A 

Escort  3 

Medium  Low 

1 

20 

A 

Escort  4 

Medium  Low 

2 

80 

A 

SAG  1 

Very  High 

3 

100 

A 

SAG  2 

Medium  high 

2 

80 

A 

SAG  3 

Medium  High 

3 

80 

B 

TBMD  1 

High 

2 

100 

B 

TBMD  2 

High 

2 

100 

B 

TBMD  3 

High 

2 

100 

B 

Escort  1 

Medium  High 

3 

80 

B 

Escort  2 

Medium  High 

3 

80 

B 

Escort  3 

Medium  High 

3 

80 

B 

Escort  4 

Medium  High 

2 

50 

B 

Escort  5 

Medium  High 

2 

50 

B 

Escort  6 

Medium  High 

2 

50 

B 

SAG  1 

Very  High 

3 

100 

B 

SAG  2 

Medium  High 

3 

100 

B 

SAG  3 

Medium  High 

3 

100 

B 

SAG  4 

Very  High 

3 

50 

B 

SAG  5 

Medium  High 

2 

50 

B 

SAG  6 

Medium  High 

2 

50 
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Table  4.  Notional  Missile  Inventory.  Nine  missile  types  were  used  for  this 
research.  The  tabled  inventory  shows  the  starting  inventory  used  for 
the  research  scenarios.  Source:  Wiederholt  (2015). 


Missile 

ESSM 

SM2 

MR 

SM2 

ER 

SM3 

SM6 

TLAM1 

TLAM2 

TLAM3 

ASROC 

Inventory 

1600 

1900 

700 

500 

400 

500 

400 

400 

1800 

Table  5.  Warplan-Mission-Missile  Requirements.  These  requirements  include 
desired  number,  minimum  number,  and  shortage  penalty  per  missile. 
Source:  Wiederholt  (2015). 


Warplan 

Mission 

Missile 

Desired  Number 

Minimum 

Number 

Missile  Shortage 
Penalty 

A 

TBMD 

SM2  ER 

40 

30 

500 

A 

TBMD 

SM3 

40 

8 

1000 

A 

STRIKE  1 

ESSM 

8 

8 

1000 

B 

TBMD  1 

SM2  ER 

40 

30 

500 

B 

TBMD  1 

SM3 

40 

10 

1000 

B 

Escort  1 

ASROC 

30 

25 

800 

B 

Escort  2 

ESSM 

16 

16 

1000 

B 

Escort  2 

SM3 

10 

0 

1000 

B 

Escort  2 

SM6 

20 

10 

1000 

B 

Escort  2 

ASROC 

30 

25 

1000 

B 

Escort  1 

ASROC 

30 

25 

800 

B 

SAG  6 

SM2  ER 

40 

20 

1000 

B 

SAG  6 

SM2  MR 

40 

40 

1000 

B 

SAG  6 

ESSM 

8 

8 

1000 
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Table  6.  Warplan-Mission  Alternate  Missiles.  The  table  lists  the  primary 
desired  missile,  alternate  missile,  and  penalty  per  missile  of  less 
desirable  type.  Source:  Wiederholt  (2015). 


Warplan 

Mission 

Desired 

Missile 

Alternative 

Missile 

Alternative 
Missile  Penalty 

B 

Escort  1 

SM2  ER 

SM2  MR 

5 

B 

Escort  1 

SM6 

SM2  MR 

8 

B 

Escort  1 

SM6 

SM2  ER 

7 

B 

Escort  2 

SM2  ER 

SM2  MR 

5 

B 

Escort  2 

SM6 

SM2  MR 

8 

B 

Escort  2 

SM6 

SM2  ER 

7 

B 

Escort  3 

SM2  ER 

SM2  MR 

5 

B 

Escort  3 

SM6 

SM2  MR 

8 

B 

Escort  3 

SM6 

SM2  ER 

7 

B 

Escort  4 

SM2  ER 

SM2  MR 

5 

B 

Escort  4 

SM6 

SM2  MR 

8 

B 

Escort  4 

SM6 

SM2  ER 

7 

B 

Escort  5 

SM2  ER 

SM2  MR 

5 

B 

Escort  5 

SM6 

SM2  MR 

8 

B 

Escort  5 

SM6 

SM2  ER 

7 

B 

Escort  6 

SM2  ER 

SM2  MR 

5 

B 

Escort  6 

SM6 

SM2  MR 

8 

B 

Escort  6 

SM6 

SM2  ER 

7 

A 

Escort  1 

SM2  ER 

SM2  MR 

5 

A 

Escort  1 

SM6 

SM2  MR 

9 

A 

Escort  1 

SM6 

SM2  ER 

7 

A 

Escort  2 

SM2  ER 

SM2  MR 

5 

A 

Escort  2 

SM6 

SM2  MR 

9 

A 

Escort  2 

SM6 

SM2  ER 

7 

A 

Escort  3 

SM2  ER 

SM2  MR 

5 

A 

Escort  3 

SM6 

SM2  MR 

9 

A 

Escort  3 

SM6 

SM2  ER 

7 

D.  COMPUTATIONAL  RESULTS 

The  three  scenarios  are  analyzed  by  observing  the  number  of  extra  missions  that 
can  be  accomplished  by  letting  the  missile  planning  tool  work  at  optimizing  the  loadout 
and  mission  assignment.  Notice  improvement  from  the  previous  case  as  we  move  from 
most  to  least  restricted  in  missile  movement  capabilities.  The  planner  advises  ship  to 
mission  assignments,  respective  ship  missile  loadouts,  and  missions  able  to  be  covered 
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given  the  restrictions  and  penalties  set  forth  in  the  model.  Each  scenario  analysis  will 
have  a  table  of  available  ships  and  each  ship’s  entry  loadout  followed  by  a  table  of 
recommended  output  ship-to-mission  assignments  and  a  table  of  missions  unable  to  be 
covered  by  the  constraints  of  the  scenario.  Heuristic  solution  quality  is  calculated  using 
the  same  relative  mixed  integer  program  gap  tolerance  calculation  as  CPLEX 
optimization  software  (2015).  This  formula  is  lbestnode-bestintegerl/(le-10+lbestintegerl) 
where  the  best  node  is  the  value  obtained  from  the  Wiederholt  (2015)  formulation  and  the 
best  integer  is  the  solution  obtained  by  the  genetic  algorithm.  Negative  integers  are 
allowed  as  input  and  the  absolute  value  operation  ensures  the  resulting  gap  is  positive. 
The  small  number,  le-10,  in  the  denominator  prevents  division  by  zero. 

1.  Scenario  I:  Fixed  Missile  Loadout 

The  first  scenario  shows  the  missile  planning  tool’s  flexibility  to  let  the 
operational  planner  keep  ship  loadouts  fixed.  This  may  be  the  case  if  missile  movements 
are  not  an  option  due  to  time  constraints,  if  the  financial  cost  of  movements  is  too  great, 
or  if  the  hazard  of  handling  is  too  high.  The  tool  will  advise  the  operational  planner  as  to 
which  missions  each  ship  should  be  assigned.  Table  7  lists  the  ships  available  in  this 
scenario  along  with  their  current  loadout.  This  scenario  includes  FDNF  ships  available  in 
cycles  1  and  2  and  deploying  ships  available  in  cycle  1.  The  value  one  in  the  “Fixed 
Foadout”  column  indicates  that  a  ship’s  loadout  is  fixed. 

a.  Scenario  I  Analysis 

After  the  overhead  time  of  building  the  sets  once  the  data  is  input  by  the  planner, 
a  solution  with  roughly  a  40  percent  optimality  gap  can  be  generated  by  the  heuristic  in 
about  10  minutes.  More  generations,  meaning  longer  run  times,  yield  between  17-22 
percent  relative  optimality  gaps.  Table  7  lists  the  fixed  loadout  of  each  ship  available.  Not 
every  ship  is  available  in  every  cycle  or  compatible  with  every  mission.  Table  8  lists  a 
resulting  ship-to-mission  assignment  recommendation  and  Table  9  lists  missions  unable 
to  be  covered  by  the  VFS  ships  as  currently  loaded.  Using  the  heuristic,  roughly  58 
percent  of  missions  remained  uncovered  during  the  fixed  loadout  scenario.  While  the 
number  of  missions  covered  in  Scenario  I  are  generally  greater  than  the  number  of 
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missions  covered  for  the  same  scenario  in  the  optimal  Wiederholt  (2015)  solution  using 
the  same  data,  the  trade-off  is  a  risk  of  imbalanced  loadouts  across  deployers. 


Table  7.  Scenario  I  Missile  Loadout.  Planners  mark  a  one  in  the  “Fixed 
Loadout”  column  to  indicate  no  desired  missile  movements  on  the 
corresponding  ship.  Each  ship  is  listed  with  its  initial  missile  loadout. 
Source:  Wiederholt  (2015). 


Ship 

VLS 

Cells 

Fixed 

Loadout 

ESSM 

SM2 

MR 

SM2 

ER 

SM3 

SM6 

TLAM1 

TLAM2 

TLAM3 

ASROC 

CG54 

122 

1 

32 

22 

18 

12 

10 

8 

8 

20 

16 

CG67 

122 

1 

32 

0 

8 

0 

1 

60 

37 

0 

8 

DDG54 

96 

1 

48 

22 

12 

3 

12 

5 

5 

5 

20 

DDG62 

96 

1 

24 

22 

12 

6 

10 

10 

0 

0 

30 

DDG56 

96 

1 

40 

2 

6 

15 

27 

0 

0 

0 

36 

DDG82 

96 

1 

40 

0 

31 

32 

0 

0 

0 

0 

23 

DDG85 

96 

1 

40 

0 

31 

32 

0 

0 

0 

0 

23 

DDG89 

96 

1 

40 

0 

31 

32 

0 

0 

0 

0 

23 

DDG60 

96 

1 

24 

20 

2 

10 

20 

0 

0 

0 

38 

DDG70 

96 

1 

48 

13 

6 

5 

14 

3 

3 

4 

36 

DDG91 

96 

1 

32 

28 

22 

20 

4 

3 

3 

4 

4 

CG65 

122 

1 

32 

30 

20 

30 

20 

6 

0 

0 

8 

CG70 

122 

1 

40 

32 

16 

10 

16 

10 

6 

6 

16 

DDG86 

96 

1 

40 

30 

20 

6 

12 

6 

4 

2 

6 

DDG92 

96 

1 

36 

21 

18 

20 

8 

4 

4 

4 

8 

38 


Table  8.  Scenario  I  Ship-to-Mission  Recommendations.  The  heuristic  output 
gives  recommended  ship  assignments  for  each  mission  in  the  fixed 
loadout  scenario. 


Warplan 

Cycle 

Mission 

Ships 

A 

cycle  1 

TBMD 

CG54,  DDG91 

A 

cycle2 

TBMD 

DDG62,  CG54 

A 

cycle2 

Strike3 

CG67 

A 

cycle  1 

Escort4 

DDG89,  DDG56 

A 

cycle2 

Escort4 

DDG82,  DDG62 

A 

cycle  1 

SAG1 

DDG86,  CG67,  CG70 

A 

cycle2 

SAG3 

DDG54,  DDG62,  CG67 

B 

cycle  1 

TBMD 

DDG85,  CG54 

B 

cycle2 

TBMD 

DDG89,  DDG56 

B 

cycle  1 

TBMD1 

DDG92,  DDG91 

B 

cycle2 

TBMD1 

DDG89,  CG54 

B 

cycle  1 

TBMD2 

DDG86,  DDG92 

B 

cycle2 

TBMD2 

DDG89,  CG54 

B 

cycle  1 

Escort  1 

DDG82,  DDG54,  DDG62 

B 

cycle2 

Escort  1 

DDG82,  CG54,  DDG62 

B 

cycle  1 

Escort2 

DDG82,  CG62,  DDG89 

B 

cycle2 

Escort2 

DDG89,  DDG82,  CG54 

B 

cycle  1 

Escort3 

DDG62,  DDG54,  DDG56 

B 

cycle2 

Escort3 

DD89,  DDG62,  DDG82 

B 

cycle  1 

Escort4 

DDG54,  DDG70 

B 

cycle2 

Escort4 

DDG89,  DDG56 

B 

cycle  1 

Escort5 

DDG60,  DDG91 

B 

cycle2 

Escort5 

DDG85,  DDG56 

B 

cycle  1 

Escort6 

DDG82,  DDG70 

B 

cycle2 

Escort6 

DDG89,  DDG62 

B 

cycle  1 

SAG1 

DDG92,  CG67,  DDG86 

B 

cycle  1 

SAG2 

CG65,  CG67,  DDG70 

B 

cycle2 

SAG2 

DDG62,  CG67,  DDG70 

B 

cycle  1 

SAG3 

CG70,  DDG62,  CG67 

B 

cycle2 

SAG3 

CG67,  DDG62,  CG54 
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Table  9.  Scenario  I  Shortfalls.  These  are  missions  unable  to  be  covered  by  fixed 
ship  loadouts.  There  are  15  in  warplan  A  and  7  in  warplan  B. 


Warplan 

Cycle 

Mission 

A 

cycle  1 

Strike  1 

A 

cycle  1 

Strike2 

A 

cycle  1 

Strike3 

A 

cycle2 

Sttikel 

A 

cycle2 

Strike2 

A 

cycle  1 

Escort  1 

A 

cycle2 

Escort  1 

A 

cycle  1 

Escort2 

A 

cycle2 

Escort2 

A 

cycle  1 

Escort3 

A 

cycle2 

Escort3 

A 

cycle2 

SAG1 

A 

cycle  1 

SAG2 

A 

cycle2 

SAG2 

A 

cycle  1 

SAG3 

B 

cycle2 

SAG1 

B 

cycle  1 

SAG4 

B 

cycle2 

SAG4 

B 

cycle  1 

SAG5 

B 

cycle2 

SAG5 

B 

cycle  1 

SAG6 

B 

cycle2 

SAG6 

2.  Scenario  II:  Fixed  and  Flexible  Missile  Loadout 

Scenario  II  adds  eight  VLS  ships  to  cycle  2  of  Scenario  I.  The  same  ships  in  the 
previous  scenario  will  maintain  their  fixed  loadout;  however,  the  eight  additional  ships 
will  have  the  opportunity  for  missile  loadout  recommendations  in  preparation  for 
deployment.  Additionally,  we  can  see  the  best  warplan-mission  assignments  for  the 
whole  complement  of  deployed  ships.  The  available  ships  and  current  missile  loadouts 
are  displayed  in  TablelO.  Notice  the  “Fixed  Loadout”  column  indicating  whether  a  ship’s 
loadout  is  fixed  or  can  be  changed. 

a.  Scenario  II  Analysis 

After  the  overhead  time  of  building  the  sets  once  the  data  is  input  by  the  planner, 


a  solution  with  roughly  an  83  percent  optimality  gap  is  available  in  about  ten  minutes. 
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Longer  run  times  have  yielded  results  as  low  as  an  81  percent  gap.  Table  10  lists  the 
loadout  of  each  ship  available.  Not  every  ship  is  available  in  every  cycle  or  compatible 
with  every  mission.  Table  11  lists  a  sample  ship-to-mission  assignment  recommendation 
and  Table  12  lists  missions  unable  to  be  covered  by  the  VLS  ships.  Roughly  54  percent 
of  missions  remained  uncovered  during  the  mix  of  fixed  and  flexible  loadout  scenario. 


Table  10.  Scenario  II  Missile  Loadout.  A  zero  in  the  “Fixed  Loadout”  column 
distinguishes  eight  of  the  West  Coast  ships  scheduled  to  deploy  in 
cycle  2.  The  loadout  planner  has  the  flexibility  to  specify  the  VLS  load 
for  these  ships  while  the  others  are  still  bound  by  the  current  loadout. 
Source:  Wiederholt  (2015). 


Ship 

VLS 

Cells 

Fixed 

Loadout 

ESSM 

SM2  MR 

SM2  ER 

SM3 

SM6 

TLAM1 

TLAM2 

TLAM3 

ASROC 

DDG1000 

80 

0 

24 

18 

7 

0 

10 

0 

21 

0 

18 

CG54 

122 

1 

32 

22 

18 

12 

10 

8 

8 

20 

16 

CG67 

122 

1 

32 

0 

8 

0 

1 

60 

37 

0 

8 

DDG54 

96 

1 

48 

22 

12 

3 

12 

5 

5 

5 

20 

DDG62 

96 

1 

24 

22 

12 

6 

10 

10 

0 

0 

30 

DDG56 

96 

1 

40 

2 

6 

15 

27 

0 

0 

0 

36 

DDG82 

96 

1 

40 

0 

31 

32 

0 

0 

0 

0 

23 

DDG85 

96 

1 

40 

0 

31 

32 

0 

0 

0 

0 

23 

DDG89 

96 

1 

40 

0 

31 

32 

0 

0 

0 

0 

23 

DDG60 

96 

1 

24 

20 

2 

10 

20 

0 

0 

0 

38 

DDG70 

96 

1 

48 

13 

6 

5 

14 

3 

3 

4 

36 

DDG91 

96 

1 

32 

28 

22 

20 

4 

3 

3 

4 

4 

CG65 

122 

1 

32 

30 

20 

30 

20 

6 

0 

0 

8 

CG70 

122 

1 

40 

32 

16 

10 

16 

10 

6 

6 

16 

DDG86 

96 

1 

40 

30 

20 

6 

12 

6 

4 

2 

6 

DDG92 

96 

1 

36 

21 

18 

20 

8 

4 

4 

4 

8 

DDG77 

96 

0 

16 

23 

13 

0 

10 

5 

15 

0 

26 

DDG90 

96 

0 

32 

20 

13 

0 

6 

23 

15 

10 

1 

DDG76 

96 

0 

32 

28 

18 

0 

9 

5 

12 

0 

16 

DDG93 

96 

0 

40 

1 

30 

4 

13 

0 

0 

0 

38 

CG73 

122 

0 

24 

18 

8 

0 

6 

12 

22 

50 

0 

DDG59 

96 

0 

24 

7 

0 

0 

3 

50 

0 

0 

30 

DDG69 

96 

0 

24 

8 

30 

0 

12 

4 

5 

0 

31 
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Table  11.  Scenario  II  Ship-to-Mission  Recommendations.  The  heuristic  output 
gives  recommended  ship  assignments  for  each  mission  in  the  fixed  and 
flexible  missile  loadout  scenario. 


Warplan 

Cycle 

Mission 

Ships 

A 

cycle  1 

TBMD 

CG54,  CG65 

A 

cycle2 

TBMD 

DDG85,  DDG69 

A 

cycle  1 

Strike3 

CG67 

A 

cycle2 

Strike  1 

DDG59 

A 

cycle2 

Escort2 

DDG76 

A 

cycle  1 

Escort4 

DDG60,  DDG82 

A 

cycle2 

SAG2 

CG73,  DDG56 

A 

cycle2 

SAG3 

DDG77,  CG73,  CG54 

B 

cycle  1 

TBMD 

DDG89,  DDG70 

B 

cycle2 

TBMD 

DDG82,  DDG54 

B 

cycle  1 

TBMD1 

DDG62,  DDG85 

B 

cycle  1 

TBMD2 

DDG85,  DDG60 

B 

cycle2 

TBMD2 

DDG89,  DDG54 

B 

cycle  1 

Escort  1 

DDG86,  CG54,  DDG82 

B 

cycle2 

Escort  1 

DDG82,  CG54,  DDG93 

B 

cycle  1 

Escort2 

DDG82,  CG65,  CG54 

B 

cycle2 

Escort2 

CG73,  DDG76,  DDG56 

B 

cycle  1 

Escort3 

DDG91,  DDG60,  DDG85 

B 

cycle2 

Escort3 

DDG93,  DDG77,  DDG90 

B 

cycle2 

Escort4 

DDG93,  DDG76 

B 

cycle  1 

Escort5 

DDG62,  DDG92 

B 

cycle2 

Escort5 

DDG1000,  DDG76 

B 

cycle  1 

Escort6 

DDG92,  DDG56 

B 

cycle2 

Escort6 

DDG69,  DDG62 

B 

cycle2 

SAG1 

DDG69,  DDG1000,  CG73 

B 

cycle2 

SAG2 

DDG69,  DDG77,  DDG85 

B 

cycle2 

SAG3 

DDG90,  CG67,  DDG69 

B 

cycle2 

SAG4 

DDG77,  DDG59,  DDG54 
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Table  12.  Scenario  II  Shortfalls.  These  are  missions  unable  to  be  covered  by  mix 
fixed  and  flexible  ship  loadouts.  There  are  14  in  warplan  A  and  10  in 
warplan  B. 


Warplan 

Cycle 

Mission 

A 

cycle  1 

Strike  1 

A 

cycle  1 

Strike2 

A 

cycle2 

Strike2 

A 

cycle2 

Strike3 

A 

cycle  1 

Escort  1 

A 

cycle2 

Escort  1 

A 

cycle  1 

Escort2 

A 

cycle  1 

Escort3 

A 

cycle2 

Escort3 

A 

cycle2 

Escort4 

A 

cycle  1 

SAG1 

A 

cycle2 

SAG1 

A 

cycle  1 

SAG2 

A 

cycle  1 

SAG3 

B 

cycle2 

TBMD1 

B 

cycle  1 

Escort4 

B 

cycle  1 

SAG1 

B 

cycle  1 

SAG2 

B 

cycle  1 

SAG3 

B 

cycle  1 

SAG4 

B 

cycle  1 

SAG5 

B 

cycle2 

SAG5 

B 

cycle  1 

SAG6 

B 

cycle2 

SAG6 

3.  Scenario  III:  Flexible  Missile  Loadout  for  All  Ships  in  All  Cycles 

Scenario  III  allows  the  planning  tool  to  adjust  each  ship-missile  loadout  in  order 
to  best  fill  mission  requirements  with  minimal  missile  movement.  This  would  advise 
operational  planners  and  Fleet  Commanders  regarding  future  deployments.  Figure  11 
shows  an  example  solved  ship  loadout.  Table  13  lists  the  final  missile  loadout 
recommendations.  Table  14  lists  the  ship-to-mission  assignment  recommendations,  and 
Table  15  shows  the  missions  unable  to  be  covered  in  this  scenario. 
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a.  Scenario  III  Analysis 

After  the  overhead  time  of  building  the  sets  once  the  data  is  input  by  the  planner, 
a  solution  with  roughly  an  88  percent  optimality  gap  is  available  in  about  10  minutes. 
Longer  run  times  have  yielded  results  as  low  as  a  78  percent  gap.  Roughly  33  percent  of 
missions  remained  uncovered  during  the  flexible  loadout  scenario  as  shown  in  Table  15. 
This  is  about  a  23-percent  decrease  of  uncovered  missions  from  the  fixed  loadouts  of 
Scenario  I. 


Table  13.  Scenario  III  Missile  Loadout.  All  zeroes  in  the  “Fixed  Loadout” 
column  allow  the  planning  tool  complete  flexibility  for  assigning  each 
ship  its  missile  loadout.  Adapted  from  Wiederholt  (2015). 


Ship 

VLS 

Cells 

Fixed 

Loadout 

ESSM 

SM2 

MR 

SM2 

ER 

SM3 

SM6 

TLAM1 

TLAM2 

TLAM3 

ASROC 

DDG1000 

80 

0 

24 

0 

28 

0 

10 

20 

0 

0 

16 

CG54 

122 

0 

32 

22 

18 

12 

10 

8 

8 

20 

16 

CG67 

122 

0 

32 

40 

8 

0 

1 

30 

23 

0 

12 

DDG54 

96 

0 

16 

40 

14 

0 

0 

30 

0 

0 

8 

DDG62 

96 

0 

24 

40 

12 

0 

0 

30 

0 

0 

8 

DDG56 

96 

0 

32 

2 

6 

15 

27 

0 

0 

0 

38 

DDG82 

96 

0 

16 

0 

40 

40 

4 

0 

0 

0 

8 

DDG85 

96 

0 

24 

2 

40 

0 

10 

30 

0 

0 

8 

DDG89 

96 

0 

40 

0 

31 

2 

0 

0 

30 

0 

23 

DDG60 

96 

0 

24 

16 

0 

10 

20 

0 

0 

0 
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DDG70 

96 

0 

48 

13 

6 

0 

14 

5 

10 

0 

36 

DDG91 

96 

0 

24 

16 

0 

20 

20 

0 

0 

0 

34 

CG65 

122 

0 

24 

38 

20 

20 

20 

6 

0 

0 

12 

CG70 

122 

0 

24 

0 

0 

0 

0 

0 

0 

0 

116 

DDG86 

96 

0 

24 

23 

22 

28 

5 

0 

0 

0 

12 

DDG92 

96 

0 

8 

21 

22 

13 

8 

0 

22 

0 

8 

DDG77 

96 

0 

24 

20 

8 

0 

16 

10 

2 

0 

34 

DDG90 

96 

0 

32 

21 

26 

0 

10 

5 

5 

5 

16 

DDG76 

96 

0 

28 

20 

12 

10 

12 

5 

5 

5 

20 

DDG93 

96 

0 

56 

20 

17 

4 

2 

4 

4 

3 

28 

CG73 

122 

0 

16 

22 

32 

3 

12 

0 

28 

5 

16 

DDG59 

96 

0 

40 

40 

0 

0 

0 

0 

30 

0 

16 

DDG69 

96 

0 

8 

0 

40 

35 

15 

0 

0 

0 

4 
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Ship  Hull 

VLS 

Cells 

VLS 

Load 

ESSM  SM2 MR 

SM2 ER 

SM3 

SM6 

TLAM1 

TLAM2  TLAM3  ASROC  Available Cells 

CG54 

122 

initial  loadout 

32 

22 

18 

12 

10 

8 

8 

20 

16 

CHANGE 

0 

0 

0 

0 

0 

0 

-8 

-9 

17 

LOAD 

32 

22 

18 

12 

10 

8 

0 

11 

33 

RISK MISSILES 

16 

0 

0 

0 

0 

0 

0 

0 

8 

DDG92 

96 

initial  loadout 

36 

21 

18 

20 

8 

4 

4 

4 

8 

CHANGE 

0 

0 

0 

0 

0 

0 

0 

0 

0 

LOAD 

36 

21 

18 

20 

8 

4 

4 

4 

8 

RISK MISSILES 

8 

0 

0 

0 

0 

0 

0 

0 

4 

DDG62 

96 

initial  loadout 

24 

22 

12 

6 

10 

10 

0 

0 

30 

CHANGE 

8 

-2 

0 

0 

0 

0 

0 

0 

0 

LOAD 

32 

20 

12 

6 

10 

10 

0 

0 

30 

RISKJVIISSILES 

16 

0 

0 

0 

0 

0 

0 

0 

8 

Figure  11.  Example  Solved  Loadout  for  Three  Ships.  Negative  values  indicate 
missile  offloads  and  positive  values  indicate  missile  onloads.  DDG92 
requires  no  change  to  its  current  loadout  to  complete  recommended 
tasking. 
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Table  14.  Scenario  III  Ship-to-Mission  Recommendations.  The  heuristic  output 
gives  recommended  ship  assignments  for  each  mission  in  the  flexible 
loadout  scenario. 


Warplan 

Cycle 

Mission 

Ships 

A 

cycle  1 

TBMD 

CG54,  DDG92 

A 

cycle2 

TBMD 

DDG62,  DDG90 

A 

cycle  1 

Escort  1 

DDG86 

A 

cycle2 

Escort  1 

DDG56 

A 

cycle  1 

Escort4 

DDG89 

A 

cycle2 

Escort4 

DDG70 

A 

cycle2 

SAG1 

DDG85,  CG67,  DDG93 

A 

cycle  1 

SAG3 

CG67,  CG70,  DDG70 

A 

cycle2 

SAG3 

DDG85,  CG67,  DDG82 

B 

cycle  1 

TBMD 

CG54,  CG65 

B 

cycle2 

TBMD 

CG67,  DDG59 

B 

cycle  1 

TBMD1 

DDG82,  CG65 

B 

cycle2 

TBMD1 

DDG76,  DDG59 

B 

cycle  1 

TBMD2 

DDG54,  DDG92 

B 

cycle2 

TBMD2 

DDG69,  DDG89 

B 

cycle  1 

Escort  1 

CG65,  DDG62,  DDG92 

B 

cycle2 

Escort  1 

DDG62,  CG54,  DDG85 

B 

cycle  1 

Escort2 

DDG60,  DDG62,  DDG92 

B 

cycle2 

Escort2 

DDG54,  DDG59,  DDG62 

B 

cycle  1 

Escort3 

DDG54,  DDG60,  DDG92 

B 

cycle2 

Escort3 

DDG65,  DDG1000,  DDG89 

B 

cycle  1 

Escort4 

DDG86,  DDG91 

B 

cycle2 

Escort4 

DDG77,  DDG69 

B 

cycle  1 

Escort5 

DDG85,  DDG56 

B 

cycle2 

Escort5 

DDG90,  DDG54 

B 

cycle  1 

Escort6 

DDG70,  DDG56 

B 

cycle2 

Escort6 

DDG89,  DDG69 

B 

cycle  1 

SAG1 

CG67,  CG70,  CG54 

B 

cycle2 

SAG1 

DDG77,  DDG76,  DDG93 

B 

cycle  1 

SAG2 

CG70,  CG54,  CG67 

B 

cycle2 

SAG2 

DDG77,  DDG62,  DDG85 

B 

cycle  1 

SAG3 

DDG92,  DDG54,  CG70 

B 

cycle2 

SAG3 

DDG1000,  DDG62,  DDG93 

B 

cycle  1 

SAG4 

DDG62,  DDG67,  DDG85 

B 

cycle2 

SAG4 

DDG90,  CG73,  DDG59 
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Table  15.  Scenario  III  Shortfalls.  These  are  missions  unable  to  be  covered  with 
flexible  loadouts  for  all  VLS  ships  in  all  deployment  cycles.  There  are 
13  missions  in  warplan  A  and  4  in  warplan  B  that  remain  uncovered. 


Warplan 

Cycle 

Mission 

A 

cycle  1 

Strike  1 

A 

cycle  1 

Strike2 

A 

cycle  1 

Strike3 

A 

cycle2 

Strike  1 

A 

cycle2 

Strike2 

A 

cycle2 

Strike3 

A 

cycle  1 

Escort2 

A 

cycle2 

Escort2 

A 

cycle  1 

Escort3 

A 

cycle2 

Escort3 

A 

cycle  1 

SAG1 

A 

cycle  1 

SAG2 

A 

cycle2 

SAG2 

B 

cycle  1 

SAG5 

B 

cycle2 

SAG5 

B 

cycle  1 

SAG6 

B 

cycle2 

SAG6 
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VI.  CONCLUSIONS  AND  FUTURE  WORK 


A.  SUMMARY 

We  set  out  to  mimic  the  data  and  penalties  of  the  Wiederholt  (2015)  formulation 
exactly.  In  this  thesis,  we  show  that  a  genetic  algorithm  heuristic  for  the  missile-to-ship 
and  ship-to-mission  assignment  problem  can  provide  a  good  and  fast  solution  for 
operational  planners.  However,  we  discovered  some  limitations  when  compared  to 
Wiederholt’s  integer  linear  program  solution  (2015).  The  greedy  heuristic  chosen  within 
the  genetic  algorithm  to  recommend  the  ship  loadouts  does  not  balance  the  missile  loads 
well  across  ships  in  the  case  that  there  is  a  large  gap  between  the  minimum  and  desired 
number  of  missiles  for  a  mission.  The  implemented  loading  heuristic  puts  the  emphasis 
on  completing  more  missions  by  not  being  short  missiles  of  the  minimum  missiles 
required  to  accomplish  a  mission.  This  puts  an  even  distribution  of  missiles  across  ships 
assigned  to  a  mission  as  a  secondary  consideration  to  mission  completion.  The 
implemented  greedy  heuristic  also  does  not  accommodate  alternative  missiles.  Future 
work  will  be  required  to  address  these  limitations.  The  Missile  Loadout  Planner  reduces 
the  uncertainty  on  how  well  we  can  be  ready  to  cover  missions  arising  from  multiple 
alternative  warplans  in  multiple  deployment  cycles  to  include  specific  mission 
requirements  and  extra  defensive  missiles  for  high-risk  missions.  Early  planning  can  also 
save  time  and  money  in  the  actual  loading  and  unloading  of  missiles  onboard  VLS  ships. 
With  the  MLP,  we  can  reduce  the  number  of  required  missile  movements  for  this  costly 
process.  The  bonus  is  that  this  tool  is  available  to  U.S.  Navy  operational  planners  on 
current  NMCI  systems  and  with  no  extra  expense  of  commercial  solvers. 

B.  RECOMMENDATIONS 

The  Missile  Loadout  Planner  is  a  flexible  planning  tool  that  allows  the  operational 
planner  to  maintain  control  of  desired  attributes  of  the  loadout  process.  The  loadout 
planning  tool’s  effectiveness  has  been  demonstrated  by  testing  against  three  scenarios 
where  the  optimal  solution  is  known.  While  the  model  was  tested  based  on  a  Commander 
Seventh  Fleet  (C7F)  Area  of  Responsibility  (AOR),  the  tool  is  flexible  enough  to 
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accommodate  any  Fleet.  This  tool,  if  implemented  for  future  ship  deployments  and 
loading  evolutions,  will  help  advise  planners  of  potential  courses  of  action.  This  will  best 
ensure  that  our  warships  have  the  appropriate  loadout  regardless  of  the  warplan-mission 
scenarios  we  face. 

C.  FUTURE  WORK 

Although  the  genetic  algorithm  heuristic  for  the  missile-to-ship  and  ship-to- 
mission  recommendations  is  a  workable  option,  more  can  be  done  to  improve  the 
effectiveness  of  the  planning  tool.  The  Missile  Loadout  Planner  could  benefit  from  being 
tested  on  real-world  data,  being  integrated  into  an  Excel  (2010)  interface,  and  refinement 
of  the  genetic  algorithm. 

1.  Real-World  Data 

The  MLP  remains  untested  on  actual  mission  requirements.  The  scenarios  in  this 
thesis  are  realistic,  but  are  based  on  fictitious  data.  We  may  discover  new  requirements 
and  restrictions  that  could  be  built  into  the  model  based  on  real-world  warplans. 

2.  Integration  into  Excel  Interface 

Work  is  already  being  done  as  part  of  the  Maritime  Operational  Planner  process 
to  tie  the  GAMS  solution  of  the  Wiederholt  (2015)  formulation  as  well  as  a  Python 
(2012)  implementation  using  the  Pyomo  package  to  a  Microsoft  Excel  (2010)  interface. 
By  adding  the  Excel  VBA  (2010)  solve  of  this  problem  to  the  interface,  planners  would 
have  three  options  to  get  loadout  recommendations. 

3.  Refine  Genetic  Algorithm  Heuristic 

While  the  Missile  Loadout  Planner  does  provide  good  loadout  and  assignment 
recommendations,  there  is  still  room  to  improve  the  heuristic.  The  loading  heuristic  used 
is  not  adequate  to  balance  loads.  More  degrees  of  freedom  could  be  obtained  by 
expanding  the  chromosome  to  incorporate  more  decision  variables.  There  is  also  room 
for  improvement  on  the  methods  for  mutations  and  natural  selection. 
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APPENDIX.  A  COMPLETE  LIST  OF  WARPLAN  MISSION 

REQUIREMENTS 


Warplan 

Mission 

Missile 

Desired  Number 

Minimum 

Number 

Missile  Shortage 
Penalty 

A 

TBMD 

SM2  ER 

40 

30 

500 

A 

TBMD 

SM3 

40 

8 

1000 

B 

TBMD  1 

SM2  ER 

40 

30 

500 

B 

TBMD  1 

SM3 

40 

10 

1000 

B 

TBMD  1 

SM6 

15 

10 

500 

B 

TBMD  2 

SM2  ER 

40 

30 

500 

B 

TBMD  2 

SM3 

40 

10 

1000 

B 

TBMD  2 

SM6 

15 

10 

500 

B 

TBMD  3 

SM2  ER 

40 

30 

500 

B 

TBMD  3 

10 

1000 

B 

TBMD  3 

mm- 

15 

10 

500 

B 

Escort  1 

ESSM 

16 

16 

B 

Escort  1 

SM3 

10 

0 

B 

Escort  1 

SM6 

20 

10 

800 

B 

Escort  1 

ASROC 

30 

25 

800 

B 

Escort  2 

ESSM 

16 

16 

1000 

B 

Escort  2 

SM3 

10 

0 

1000 

B 

Escort  2 

SM6 

20 

10 

1000 

B 

Escort  2 

ASROC 

30 

25 

1000 

B 

Escort  3 

ESSM 

16 

16 

1000 

B 

Escort  3 

SM3 

10 

0 

1000 

B 

Escort  3 

SM6 

20 

10 

1000 

B 

Escort  3 

ASROC 

30 

25 

1000 

B 

Escort  4 

ESSM 

16 

16 

1000 

B 

Escort  4 

SM3 

10 

0 

1000 

B 

Escort  4 

SM6 

20 

10 

1000 

B 

Escort  4 

ASROC 

30 

25 

1000 

B 

Escort  5 

ESSM 

16 

16 

1000 

B 

Escort  5 

SM3 

10 

0 

1000 

B 

Escort  5 

SM6 

20 

10 

1000 

B 

Escort  5 

ASROC 

30 

25 

1000 

B 

Escort  6 

ESSM 

16 

16 

1000 

B 

Escort  6 

SM3 

10 

0 

1000 

B 

Escort  6 

SM6 

20 

10 

1000 

B 

Escort  6 

ASROC 

40 

25 

1000 

B 

SAG  1 

SM6 

10 

6 

600 

B 

SAG  1 

TLAM1 

30 

30 

800 

B 

SAG  1 

SM2  ER 

40 

20 

1000 
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B 

SAG  1 

SM2  MR 

40 

40 

1000 

B 

SAG  1 

ESSM 

8 

8 

1000 

B 

SAG  2 

SM6 

10 

6 

1000 

B 

SAG  2 

TLAM1 

30 

30 

1000 

B 

SAG  2 

SM2  ER 

40 

20 

1000 

B 

SAG  2 

SM2  MR 

40 

40 

1000 

B 

SAG  2 

ESSM 

8 

8 

1000 

B 

SAG  3 

SM6 

10 

6 

1000 

B 

SAG  3 

TLAM1 

30 

30 

1000 

B 

SAG  3 

SM2  ER 

40 

20 

1000 

B 

SAG  3 

SM2  MR 

40 

40 

1000 

B 

SAG  3 

ESSM 

8 

8 

1000 

B 

SAG  4 

SM6 

10 

6 

1000 

B 

SAG  4 

TLAM1 

30 

30 

1000 

B 

SAG  4 

SM2  ER 

40 

20 

1000 

B 

SAG  4 

SM2  MR 

40 

40 

1000 

B 

SAG  4 

ESSM 

8 

8 

1000 

B 

SAG  5 

SM6 

10 

6 

1000 

B 

SAG  5 

TLAM1 

30 

30 

1000 

B 

SAG  5 

SM2  ER 

40 

20 

1000 

B 

SAG  5 

SM2  MR 

40 

40 

1000 

B 

SAG  5 

ESSM 

8 

8 

1000 

B 

SAG  6 

SM6 

10 

6 

1000 

B 

SAG  6 

TLAM1 

30 

30 

1000 

B 

SAG  6 

SM2  ER 

40 

20 

1000 

B 

SAG  6 

SM2  MR 

40 

40 

1000 

B 

SAG  6 

ESSM 

8 

8 

1000 

A 

STRIKE  1 

ESSM 

8 

8 

1000 

A 

STRIKE  1 

TLAM3 

50 

50 

1000 

A 

STRIKE  2 

ESSM 

8 

8 

1000 

A 

STRIKE  2 

TLAM2 

50 

50 

1000 

A 

STRIKE  3 

ESSM 

8 

8 

1000 

A 

STRIKE  3 

TLAM1 

50 

50 

1000 

A 

Escort  1 

ESSM 

8 

8 

1000 

A 

Escort  1 

SM2  ER 

40 

30 

1000 

A 

Escort  1 

SM6 

20 

10 

1000 

A 

Escort  1 

ASROC 

30 

25 

1000 

A 

Escort  2 

ESSM 

8 

8 

1000 

A 

Escort  2 

SM2  ER 

40 

30 

1000 

A 

Escort  2 

SM6 

20 

10 

1000 

A 

Escort  2 

ASROC 

30 

25 

1000 

A 

Escort  3 

ESSM 

8 

8 

1000 

A 

Escort  3 

SM2  ER 

40 

30 

1000 

A 

Escort  3 

SM6 

20 

10 

1000 

52 


A 

Escort  3 

ASROC 

30 

25 

1000 

A 

Escort  4 

ESSM 

8 

8 

1000 

A 

Escort  4 

SM2  ER 

40 

30 

1000 

A 

Escort  4 

SM6 

20 

10 

1000 

A 

Escort  4 

ASROC 

30 

25 

1000 

A 

SAG  1 

SM6 

10 

6 

1000 

A 

SAG  1 

TLAM2 

30 

30 

1000 

A 

SAG  1 

SM2  ER 

40 

20 

1000 

A 

SAG  1 

SM2  MR 

40 

40 

1000 

A 

SAG  1 

ESSM 

8 

8 

1000 

A 

SAG  2 

SM6 

10 

6 

1000 

A 

SAG  2 

TLAM2 

30 

30 

1000 

A 

SAG  2 

SM2  ER 

40 

20 

1000 

A 

SM2  MR 

40 

40 

1000 

A 

ESSM 

8 

8 

1000 

A 

SM6 

10 

6 

1000 

A 

TLAM2 

30 

30 

1000 

A 

SM2  ER 

40 

20 

1000 

A 

SAG  3 

SM2  MR 

40 

40 

1000 

A 

SAG  3 

ESSM 

8 

8 

1000 
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